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■ Abstract. We propose three effective Hamiltonians which approximate atoms 

f - ") in very strong homogeneous magnetic fields B modelled by the Pauli Hamilton- 

f"^ ' ian, with fixed total angular momentum with respect to magnetic field axis. 

, All three Hamiltonians describe N electrons and a fixed nucleus where the 

Coulomb interaction has been replaced by independent one-dimensional effec- 
tive (vector valued) potentials but without magnetic field. Two of them are 
solvable in at least the one electron case. We briefly sketch how these Hamil- 
tonians can be used to analyse the bottom of the spectrum of such atoms. 
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1. Introduction 

The Pauli Hamiltonian of a non-relativistic atom with an infinitely heavy nucleus 
and electrons with spin in a constant magnetic field B of strength B is given by: 



in 

(N 

> 

o 

^ ■ N n i i 7 \ 1 

O. (!) i^-iBA^.B-^ + JJ ^ 

j=l v 1 Jl/ l<j<k<N 1 J B| 

where rj = (xj, yj, Zj) £ R 3 are the coordinates of the j-th electron, a 3 - is its spin, 
and Vj is the gradient with respect to rj. Note that we have made the choice of 
^B A r for the vector potential of B, and that we are working in atomic units. 

We fix the direction of B to be the z-direction: B = 5(0,0,1) with B > 
without loss of generality. Recall that the z-component of 5j is given, in the Pauli 



o 

is 

Oh 

representation, by 

1 ( -1 o 

1 • • • a Zj ■ ■ ■ ® /, a Zj = - I Q 1 



acting on the iV-fold tensor product C 2 ® • • • <E> C 2 . It is known (see [KaKuj ) that 



H B := H t! (Z, N) defines an essentially self-adjoint operator on 



N 

2 /td)3\ ,o< rr<2 



(2) H = (X)L 2 (W i )(g>C 



3- 



the Hilbert space of distinguishable electrons or "boltzons" (meaning particles sat- 
isfying the Boltzman statistics: we thank Beth Ruskai for introducing us to this 
expression) . Physical atoms are of course modeled by H B restricted to the fermionic 
subspace 

jv 

(3) Ht = /\i 2 (K 3 )(g)C 2 , 

of totally anti- symmetric wave-functions in 7i, where A stands for exterior product. 
A useful alternative description of 7i, used in atomic physics, is given by the unitary 
map 

H -> L 2 ((R 3 x {±1})*) , ^ -» U(tP)(n, si; ■ ■ ■ ; r Nl s N ), 
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the isometry being defined by taking the components of ip(ri,- ■ ■ , rjv) with respect 
to the natural basis of C 2 <£> • • • (8 C 2 consisting of the products ®je Sj of the nor- 
malized eigenvectors e± of a Zj , e_ corresponding to "spin down" and e + to "spin 
up". In this new representation, a Zj acts as the multiplication operator by Sj/2, 
and the fermionic subspace Tlf of (J2J) is simply obtained by anti-symmetrizing with 
respect to the 4-tuples of variables (r±, si), • • ■ , (rjv, sjv). 

The main results of this paper can be summarized in the following five theorems 
below. It is worthwhile to observe that H B commutes with each individual spin 
operator a Zj and therefore decomposes in a direct sum which is unitarily equivalent 
to 

e i^-- + £(i+^.)§, 

where S z :— Ylj—i & Zj denotes the z component of the total spin operator. So from 

now on we will consider only H B ' S * = ~ I V~ and denote it again by H B . Notice that 
this operator simply acts in (g)^ =1 L 2 (R 3 ). Since the Hamiltonian H B also commutes 
with the total angular momentum operator in the field direction, which we will call 
L z , we can fix a value M > of the latter. Our results will imply that the bottom 
of the spectrum of H B will necessarily occur for a non-negative value of M, and we 
will therefore restrict ourselves to M > 0. Let H B ' M be the restriction of H B to 
the M-th angular momentum channel (in the field direction), and let n2j M be the 
orthogonal projection onto the lowest Landau states with z-angular momentum M 
(cf . for the precise definition) . We define the effective Hamiltonian h c s — h^ M 
by 

(4) = nf ff M h b <™ nf ff M , 

and let 

ttB.M _ T YfB.U 
LL ± — 1 ll cff 

be the projection onto the orthogonal complement (always restricting ourselves to 
the M-th z-angular momentum channel). The operator h^g M is the first, and most 
encompassing, of three "effective Hamiltonians" we will consider in this paper. The 
two other ones, called h^' M and hf , will be defined below. It will be convenient to 
complete /i^ M as follows 

/r\ tjBM u BM^ttB,M tj-B.M n fl,M rrfl,M n S,M 

(5) H cS := hj ®H ± ' , with H ± :=U ± ' H ' J TI X ' . 

For any self-adjoint operator A, we let a(A) denote the spectrum of A and p{A) its 
resolvent set. Our first main result is: 

Theorem 1.1. Let a — a{B) be the unique positive solution of the equation 

(6) a + log a = -logB, 

and let d ff(£) = dist c(h^ M )j . There exist positive constants B c g, c c r and 

C e s, which only depend on Z, N and M, such that for all B > B c g, and all real £ 
satisfying 

a . . 1 2 

Coff— = < 4ff(£) < 2 a ■ 



a{Bf 



we have that £ G p(H B ' M ), and 

(7) ||(^,M_ -l_ (ff B,M_ crl| < Ceff 
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Remarks 1.2. (i) We shall see in Theorem IXD that a(Hf M ) C M+ for B > B^, 
where the latter is defined by formula J3HJ 1 . Since B c s > B ffift by p()7(l . and 
since one can see from © and (HHJl below that R+ C a(h^ ), it follows that 
cr(/i c g ) = cr{H cS ) when B > B eS . 

(ii) The equation for a(B) is equivalent to ae a — \/~B, and can therefore be written 
as 

a{B) = W{\fB), 

where W is the principal branch of the Lambert W function, see e.g. [CoGoHaJeKn . 
Using known properties of the Lambert W-function, or by elementary arguments, 
one shows that 

(8) a (i?) = Ilog J B-log( 2 )i? + log2 + 0^^i^^ , 5^ co, 

where log^(a;) :— log(logx), for x > 1. In particular, a(B) ~ log(v A S), as B — ► oo. 

(iii) Our proof yields explicit constants i? G ff, c e ff and C e s, for given N,M. and Z. 
This is also true for the constants in Theorems 11.31 and 11.51 below. 

(iv) The upper bound a 2 / 2 on d c s (£) is only there to allow a simple expression for 
the upper bound in J7J), and is by no means essential. The same remark applies to 
Theorems O ESI and O below. 

Some applications of Theorem 11.11 as well as of Theorems 11.31 11.51 ll.fil and 11.81 

below, to the study of spectral properties of H B ' M are given in the concluding 
remarks section. 

The operator hfg M has the structure of a multi-particle Schrodinger operator on 
the real line, R: 

(9) Ks =-2 A ~ Z Z^ V j y z i> + l^ V jk ( z j~ z k), 

3 ]<k 

with operator-valued potentials acting point-wise on a certain finite-dimensional 
Hilbert-space F^, defined in l|32|) below. Essentially, is the vector space spanned 
by the lowest Landau states with angular momentum L 2 = M . As we will see, 
Ran II^ M is canonically isomorphic to the space L 2 (M. N , F^) of vector- valued L 2 - 
functions , and the potentials in © are simply obtained by projecting the respective 
Coulomb terms in along Ib^ M : 

(10) V. ' ( Zj ):=U e 4 t—^ c s V jk fe-^) :=n efi TT—-\ n cs ■ 

I ' 3 I I ' j ' k I 

We will show that the potentials i|l(J|) can be approximated by certain simpler ones, 
which will give rise to our two other effective Hamiltonians. Define the tempered 
distribution q B on M by: 

(11) q B (z) =log5<5(z) + Pf ' ' 
where 

K (r) : = ^Csgn(^)log|^|), 

(with distributional derivative) is the finite part (in the sense of Hadamard) of the 
singular function l/|z|; Pf(|z| _1 ) should be interpreted as a regularization of the 
(3-dimensional) Coulomb potential restricted to the line. Also introduce (constant) 

^see our convention on constants at the end of this introduction 
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finite dimensional operators C™' B ^ M an d Cjjf' M , acting on the vector space 
introduced above, denned by 

,BM ._ ttB.M, ( -° f 2 | 2\ \ tt-B,: 



( 12 ) C J : =- n cff lQ g( jfcj +Vj) j n cff . 



and 

no\ ^,e,B,M tt-B, Mi (" ll \2 , l \2\ \ t-tB, 

(13) := -II eft log \ — {{xj - a*) + {y 3 ~ y k ) )Jn ef j . 

The superscripts " n" and " e" stand for " nucleus" and " electron" , respectively, as a 
reminder that l|12f) is a vestige of the interaction between the j-th electron and the 
nucleus, while (|13(l originates in the electron-electron interaction between electrons 
j and k. Finally, define an operator h B ' M on L 2 (R, F^) (the 'C standing for 
'Coulomb') by 

(14) h*» = -\a-zJ2(i b ^) + c;^5( Zj ) 

j 

+j2{i B ^-^) + c;f M 5(z J -z k ) 

j<k 

As we will see in section 4, the right hand side of 1)14(1 defines a self-adjoint operator 
h c ' on L 2 (R Ar ), despite the distributional potentials. This will be a consequence of 
the Kato-Lax-Lions-Milgram-Nelson Theorem. The form domain of h^' M is simply 
the vector- valued first Sobolev space iJ 1 (]R Ar ; Fm), while its operator domain will 
be characterized in appendix A. As in (jSJ we introduce 



tjBM ,B,M „ rr-B.R 

H c ' := ti c ' © H ± ' 
Our second main theorem then is the following: 



Theorem 1.3. Let a — a(B) be as in Theorem ] and put dc (£) := dist \ £, cr(h^' 

There exists positive constants Be, cc and Cc which depend only on Z,N and M, 
such that for all B > Bq and all real £ satisfying 

o?l 2 1 

we have that £ E p(H B ' M ), and 



(15) || (H B > m - g)" 1 ~ (gg' M ~ 0" 1 II < 



Remark 1.4. In top order in B, all that remains of the electrostatic potentials 
in H B ' M are the extremely short-range 5-potentials. In next order, the long range 
character of the original Coulomb potentials reasserts itself in two ways: in the mag- 
netic field direction, through the Pf (| • | -1 )-terms in h^' M , and in the transversal 
directions, through the (j™' B,m - and C^ B ' M - terms. The latter are in fact sim- 
ply the quantum mechanical mean, with respect to the projection onto the lowest 
Landau band states of total angular momentum M (in the field direction), of a 
2-dimensional logarithmic potential, minus a independent constant. This logarith- 
mic potential is the natural electrostatic potential for the plane. Physically, this 
can be understood as follows: under the influence of the strong magnetic field the 
electrons will spiral closely around the field lines, along circles of radius 0{B~ 1 / 2 ) in 
the plane transversal to the field, while occupying an interval of size 0((logB) _1 ) 
in the field direction itself, as a consequence of the nuclear attraction. For big 
B, and at different locations in the (x, y)-plane, they will see each other and the 
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nucleus as so many infinitely long charged wires, and, as is known from classical 
electrostatics, such wires interact via a logarithmic potential. 

A simpler effective Hamiltonian, our third and last one, and historically the first 
to be proposed (cf. |LSY| . [BaSoYj . |BD) ). is roughly speaking obtained by only 
keeping the leading term in the potential of /i^' M . More precisely, we put: 

(16) hf = -\& z +2a{B)v s , 



where 



N 



(17) v s (z)=-Zj2S(z j ) + ^2S(z j -z k ). 

3=1 j<k 

Looking back at (jlljl it would seem natural to take as potential log Bvs, but it 
turns out that 2a(B)v$ leads to smaller error estimates; notice that in view of (JSJ, 
2avg is also a 0(logi?) part of the potential in hf,. Furthermore, with this choice 
the coupling constant 2a(B) is positive for all B > which is not the case for log B. 
Contrary to our previous two effective Hamiltonians, hf does not explicitly depend 
on M anymore, but it will operate on an M and B-dependent Hilbert space, namely 
L 2 (M. N , F^) ~ L 2 (Ht N ) (gi (which in fact are canonically isomorphic for different 
B). Considering ((TBI) as acting on scalar L 2 (R N ), we define the 8- model as being 
the operator 

(18) hf' M := hf (g> IpB, I f b being the identity operator. 

We will often simply write hf, except when we want to stress the vector- valued 
nature of the L 2 -functions in the domain. Again as in Q we introduce 

H s := h 5 ' ®H X > . 
Our third approximation theorem is: 

Theorem 1.5. Let a := a(B) be as in Theorem } 1. 11 and put ds{^) '■= dist o~(hf' M ) 
There exist positive constants Bg, c$ and Cg, depending on N , Z and M, such that 
for all B > B$ and real £ satisfying 

(19) c s a < d s (0 < ia 2 , 
we have that £ S p(H B ' M ), and 

(20) || (H B <™ - 0" 1 - (H^ M - O-'W < 

See jBDj for weaker versions of this theorem. We also mention [BaSoYj . which 
established 2 the convergence of the ground state energy of fermionic H B (see be- 
low) to that of bosonic (scalar) hf on L 2 (R N ), using variational arguments: these 
authors did not fix M, but they only proved convergence of the ground state energy, 
while we can conclude much more from the norm resolvent convergence to the effec- 
tive hamiltonians; see §9 for a list of applications of the results of the present paper. 
Earlier, jLSYj had shown that the ground state of the Hartree mean- field model 
associated to (|lfit approximates the quantum mechanical ground state energy in 
the so-called hyper-strong limit Z,B/Z 3 — > oo, assuming N/Z uniformly bounded. 
The idea that a model such as the 8- model could be relevant in the context of strong 
magnetic fields is not new in the physics literature, see e.g. |Spr|. 



2 fBaSoYI . following ILTT7I . first did a re-scaling of H B 's ground state energy which allowed 
them to compare with hf for a fixed B (e.g. fixing 2a(B) = 1). Since this homogeneity property is 
not valid anymore for our other two effective hamiltonians, we prefer not to do this here (contrary 
to our earlier papers IBDI ). in order to have a coherent presentation. 
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We next turn to the effects of particle symmetry. Electrons in physical atoms are 
fermions, and we now consider the analogues of Theorems 11.11 IP1 and n~5l for H B 
restricted to the fermionic subspace Hf = P AS (H), where P AS is the orthogonal 
projection onto the subspace of anti-symmetric symmetric wave-functions defined 
by: 

P AS ip(ri,Sx;--- ;r N ,s N ) := — ^ ("^V (Mi) > s t(i) > ' ' ' > r «(N), s a (N)) , 

o-eSiv 

with (-1) CT := (-l) s § n The projection P AS commutes with H B , L 2 , § z and 
with the TV-particle Landau Hamiltonian Hq defined in l|24|) below, and there- 
fore also with II « (see section 2). Recalling that we have fixed our spins to 
§ z = —NB/2, P AS for us will only act in the 'spatial' variables (rj., • • • , r^). Let 
H B,M ,_ p AS H BM P AS = h bm p as^ the f ermionic Pami operator with z-angular 

momentum M. Similarly, introduce 'fermionized' versions of the other operators: 

tt-B,M ._ pASTT-B.M -prB.M ._ p ASttB,M ttB,M ._ p ASttB,M _ pA S >B,M o-B,M _. 
n eff,f ■- ^ n cff ' U _L,f •- ^ U _L > -"off.f •— ^ U cS - ^ n eS ® U ±,f ~ • 

ftg* © H*f, h*? = P As h B c M , hf A := P As hf, = P AS H^ M and finally 

-ffff M := P AS H^' M . A careful examination of the proofs of Theorems II . II fOl and 
II .5! will show: 

Theorem 1.6. Theorems \l.l\ M.'A and \1.5l also hold true for the fermionized oper- 
ators. For example, if ^ satisfies the conditions of theorem \l.l\ with d e ff(£) replaced 

by deff,f(0 :=dist U,(r(h^f), then 

) -1 rv(FS) 2 

*off,f 



(21) || (V' M - £ 



u. L ;,\~-a ii en- n[B) ~ 



with the same constants as before. Similarly for , Hf^. 

Remark 1.7. Theorem ll.tjl is not simply obtained by "sandwiching" Theorems ll.il 
11.31 and ll . 51 between P AS , since the statements thus obtained would not involve the 
distances to the spectra of the fermionized operators. Also, we established the 
fermionic versions with the same constants as for the boltzonic ones, but it is 
conceivable that one could have smaller constants in the fermionic case. 

The operators h^A f and h 3 ,'^ are easily described: 
1 N 

(22) ft ef j f = --A - F av ; x (Zj) + 2^ V av : 2 { Zj - z k ), 

j=l j<k 

where, using the notation of (|10fl . 

t/B,M/ \ 1 \~^ T/ B,M/ \ 

j 

and 



t/B,M/ \ I N \ >T^ T/ _B,M/ n 



i<fc 

B,S 



2 

the average one-, respectively two-particle potentials. Similarly, ftp'™ equals 
(23) ft c> f = -iA-^E(^^) + ^ M ^) 



,b,: 



- ( 9 S ( z j ~ z k) + Cwf 2 ,M 5(2!j - z k ) 
j<k 
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with 



°av:l 



and 



rf e,B,M _ ( N \ ST^ n e,BM 
°av:2 — I 2 i Z^Sfe ' 
V 7 j<k 

while hff 1 is given by the same expression as hf , except of course that its domain 
changes. 

To complete the picture, we finish with an explicit description of Ran(n^^), 
which is equal to P AS (Ran IL^ M ) = P AS (L 2 (R N ) ® Fix M > 0, and let 

S(M) = {m £ N N : \m\ = M}, 

the set of partitions of M, where, as usual, \m\ = mi + - ■ -+rriN if m = (mi, • • • , TOjy). 
The vector space is spanned by the lowest Landau (generalized) eigenstates in- 
dexed by to G S(M), cf. section 2 below. The symmetric group Sm acts on S(M) 
in the natural way, by permuting the indices of an element m = (mi, • • • , Ton) of 
S(M). Under this action, £(M) will decompose as a finite union of disjoint orbits: 

S(M) = |J S N ■ to, 

me M 

M. being a set of representatives of the orbits. If denotes the stabilizer of 
to e M, then we write M as a disjoint union M, = M\ UM2, with M.\ the subset 
ofthosern € M such that G s = {e}, and A^2 = M\Mi, its complement. In other 
words, to £ iff no two components of to are the same, and to G .M 2 iff at least 
two of its components are identical. Let L\ S (HL N ) be the space of anti-symmetrical 
wave functions in L 2 (R N ). We then will prove, in section 8, that: 

Theorem 1.8. There is a natural unitary isomorphism 



e 

U$ : Ran(lO -> L 2 (R N ) JT v ^ 



and 



Remark 1.9. The operators /i^'f U^* an< ^ f 1 ^m* can m i x different 

components of Ran (U^), as we will see at the end of section 8, and will therefore 
have a more complicated structure. 

The paper is organized as follows. Section 2 contains the precise definition of our 
effective projector n^j M . In section 3 we establish, with the help of the Feschbach 
decomposition, a first approximation theorem, comparing H B ' M, s resolvent at £ 
with that of iQ M + W S ' M (£), where the last term is an auxiliary 'potential' which 
itself depends on the spectral parameter £. Section 4 analyzes the large- B behavior 
of the potential of H^ M , as well as that of W B ' M . Sections 5, 6 and 7 are devoted 
to the proofs of, respectively, theorems 11.11 11.31 and 11.51 In section 8 we prove 
theorems 1 1 . 61 and 1 1 . 81 Section 9, finally, concludes with some applications to the 
spectral theory of H B ' M , and some general observations. 

Convention on constants. In the course of this work, we have had to introduce 
a large number of constants. To keep track of them, we will use the convention that 
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whenever the subscript of a constant is a number, the number refers to the formula 
where the constant in question was first introduced. That is, Ci x \ :— constant 
defined in formula (x). 

Acknowledgements. We first of all would like to thank Beth Ruskai for having 
introduced us to the fascinating world of atoms in strong magnetic fields. We would 
also like to thank the organizers of the Clausthal PDE Conference, Summer 1 999, 
and of the Critical Stability Workshop at les Houches, October 2001, for having 
given us the opportunity to report on preliminary written versions of this work, 
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2. Non-interacting electrons and the Lowest Landau Band 

We begin by reviewing the spectral decomposition of 

N 1 / 1 1 

(24) H* := H Q B (N) := £ - (-V„ - -B A r 3 ) 2 - NB 

3 = 1 V 

the "free" Hamiltonian of N independent electrons interacting only with the field B. 
Recall that we have fixed all electron spins in their s Zj = — 1-state. The operator 
Hq is just a direct sum of N one-particle operators | ((j _1 V r - |BA r) 2 — B), 
whose spectral decomposition is explicitly known: 



(\pI + J (2n + |m|-m))n 



B 



Here, p z is the momentum in the field direction, and n is the projection, in the 
x, y-variables, onto the normalized eigenfunctions Xmn = Xmn( x >y) G L 2 (M. 2 ) of 
the operator 

1 a B 2 9 o, B 
--A x , y + —(x 2 +y 2 )--, 

restricted to the m-eigenspace of L z = xp y — yp x , the angular momentum in the 
field-direction. These eigenfunctions are explicitly known in terms of Laguerre 
functions, see e.g. .1 W I. but for our purposes we will only need those with n = 0, 
m > 0. These have a particularly simple expression: if (p, tp) are polar coordinates 
in the x, y-plane, then 

(25) Xl := Xl, : {x,y) - ^-e^. 



The spectral decomposition of H ( f is simply the sum of the one-particle de- 
compositions, and the projections onto its eigenstates will be indexed by TV-tuples 
m = (mi, • • • j ttin) G Z , n — (n\, ■ ■ ■ un) G N^. If we let 



then 

(26) = m I > -<r +^ ( 2,v,-|m, -,„,) j n;i„. 



m,n \ j — 1 



The Lowest Landau Band of Hq is defined as 

(27) = Ran B^ , 
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where := (0, • • • ,0), and if we put 



N 



(28) X%(x,y) := Xm 4 (^j . ) . mi ' " " " ' mjv - °' 

then £ will be spanned by the tensor products X£ (g> u, with w = u(z) e L 2 (M Ar ). 
We will call the X B the lowest Landau band states (these are not eigenvectors of 
Hq , but X^®1 would be generalized eigenvectors with eigenvalue 0). The operator 
Hq restricted to C B simply is the free Laplacian in the field direction, 

2 2-^ p zj ~ 2 

3 

where z = (z\, • • • , zn). 

We next reduce the Hamiltonians H B and H B to their angular momentum 
sectors with respect to the field direction. The total orbital angular momentum 
in the direction of B = (0, 0, B), 

L z = ^(XjPyj ~ VjPxj), {Px,Py,Pz) = \^r, 
3 

commutes with H B and Hq , and is therefore a constant of motion for both Hamil- 
tonians. If P M is the orthogonal projection onto the M-th eigenspace of L z , then 
we let 

(29) H B - m := H B P m , << M := H B P u 

-2fT®3N\ ^ {n2N 



acting on H := L (M? ) <g> C . Since we are primarily interested in the spectral 
behavior of H B near the bottom of its spectrum, we will restrict M to Z + , for 
C B n Ran P M ± {0} ^ M > 0. Indeed, notice that since H B " M is unitarily 
equivalent to H B ' M + MB, one has inf a(H B |m>o) < hif <t(H b \m<o) as soon as 
B > 0. 

We next let 

(30) S(M) = {m= (m 1 ,---m N ) € Z N : m d - > 0, m 1 + ---+m N = M}, 
the set of partitions of M, and define the effective projection n^j M by 

(31) nf ff M := V U B m 



eff ' — 1L m,Q- 

m£S(M) 



This is simply the orthogonal projection onto Co n (L z = M}. We also let II ± 



BM 



be the orthogonal projection onto the orthogonal complement of Ran(II cf j ) in 
Ran(P M ). Observe that 

nf ' M n 



mi + . . . + mjv = M, 
J2j 2nj + \m,j\ — rrij > 2 

If we let be the finite dimensional vector space spanned by the lowest Landau 
states with total angular momentum M, 

(32) F B := Span{X* : m e S(M)}, 

then we can identify the range of II^ M with the space L 2 (W N , F^) of F^ -valued 
X 2 -functions, as we will do without further comment. 

To lighten the notations, we will often suppress one or both upper-indices B or 
M, unless where this would cause confusion. This will always be clearly indicated, 
usually at the beginning of a section. 
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3. Estimates for Feschbach decompositions 



We fix a non-negative integer M > 0. In this section we will drop all upper- 
indices B, M, and simply write H for H B ' M , Hq for Hq' M and II e ff respectively II j_ 
for n^ M and 1Tj j,m . We write our atomic Hamiltonian H B,M as 

H = Hn + V, 



where 

J J 3<k J 

the electrostatic potential, and introduce the operators 

v cff := n off vn cff , v± := nxvrix, v±, bS := nxvn cfr , 

and its adjoint, V e ff,x = n e ffVIIx- These are to be considered as operators on 
Ran Hgg, Ran H±, and between these two Hilbert spaces, respectively. We further- 
more put 

h e e := n c ffi7n c ff , H± := II_i_iiTI_i_, i?j_, e ff := n_i_.HTIeff = H* s ± , 
and for (eC introduce the resolvents (wherever defined) 

R := R(0 := (H ± » «3f := 00 := (^eff + W(£) - 0~ X , 

where 

W := W(0 = -Veff,xfl(OVx,eff. 

l off 



Strictly speaking is not a resolvent since the potential W(£) depends on the 
spectral parameter £. The operators R and i?^ act on, respectively, the ranges of 
ITx and of II c ff . Finally, let 

T := H P M , T eS := n off Tn cff , Tj_ := n ± TII ± ; 

T commutes with H e g and flj^, and n^TIToff = 0. Note that 

T± >Mi, 

on the range of P M . 

Using matrix notation associated to the decomposition p u H = Ran II e ff © 
Ran 11^ , we decompose H as 

h c « H c ff,± \ _ I T c ff + V c fF V c ff ,j_ 
#x,eff H± J-\ V ± , off Tf+Vx) 

By the classical Feschbach formula, we then have 



(34) H 



(35) {H-0~ l = 



-RV ±jeS R^ R + RV^ cS R™V cS ^R 



for those £ 6 C for which the right hand side makes sense. The following theorem 
is the main result of this section: recall that p(A) denotes the resolvent set of an 
operator A, and a(A) its spectrum. 

Theorem 3.1. Let 

(36) B^:=16Z 2 N(M + N + 2), 
and 

c 2 

( 37 ) C m ■= c o + i ° . c o = (32Z 2 7V + 8iV(7V _ l) 2 )(M + iV + 2). 
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If £ < an d if the field strength B > B [^fi[ , then £ € p(H±). If, moreover, 
£ £ a(H eS + W), then £ G p(#), and 



(38) 



ii (ir - er 1 - (fteii + w - er 1 © i?(on < 



c, 



E3_ 



where d^(^) is the distance of £ to o~(h e g + W). 

Proof of Theorem \S.l\ The proof consists of estimating the relevant matrix elements 
in the Feschbach formula. This will be done in several steps. Let 



J2o = (Tx-0~ 1 , 



the resolvent of Tj_ on Ran II j_ • 



Bound on i?o- Since T± > B on Ran II j_ and since £ < 0, it immediately follows 
that ||i? || < B~ x . Write V as 

(39) V = V„ + V e , where V„ := - 5") A, and V e := V ^ r, 

j l r il j<k \ r J rfc l 

the sum of electron-nucleus and the electron-electron interactions. 

A remark on notation: we will often leave the projection understood when 
multiplying operators on the left and/or right by Rq or R, and for example simply 
write Rq V n RQ instead of the more explicit Rl^IixVrJi^Rl' 2 . 

Bound on i?o /2 V„i?o /2 . First, since Rq < B on the range of II , 



< (v^RoK 



Rq) 2 < B 1 V RnV 2 \/ Rn 



By Cauchy-Schwarz, 



^V 2 ^/W < Z 2 Ny/R^Yl \ r i\~ 2 VRo 



Next, using Hardy's inequality on 



< — 4A,-, and the fact that 



H fl =^-iA j + i|BAr i | 2 -f(L, 



N), 




A,- 



/?„(rL + |(j 



TV) 



'i?,n 



= 8n_L + 8£i? + 4(M + N)BR 
< (8 + 4(M + 7V))n ± , 

since £ < 0. It follows from these estimates that 

||i^ /2 V„i^ /2 || < 2ZB-^ 2 ^N(M + N + 2). 



We note, as a consequence, that if 6 = 4Z 2 7V(M + iV + 2), then \\Rl /2 V n R^ ' \\ < 
(boB -1 ) 1 / 2 < 1 if B > bo- For later reference we also note the: 

Bound on Rq^ 2 V 2 Rq^ 2 : the estimates above immediately imply that this positive 
operator is bounded from above by 4Z 2 N(M + N + 2). 



?V 2 i 
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Rn 



Existence of and bound on R. Since the electron-electron repulsion V e > 0, it 
follows that R < i?Ni, where i?Ni = (T± + V n ,j_ — the resolvent of an atom 

with non-interacting electrons. Using the symmetrized resolvent formula, 

(40) i?Ni = VR~o (l + V^Vn^Rt 

we see that i?Ni exists and is positive if B > bo and £ < 0. Hence T± + V n .± — £ > 
and therefore also H± - £ and R. Moreover, if B > 4b = 16Z 2 N{M + N + 2) = 
B ffifi[ , then every £ < belongs to p(H±) and 

0<R j^ 11 <B~ 1 ( ,. 1 K2B- 1 



„ ^ /2 v r ^ /2 || " Vi-^os- 1 ) 1 / 2 

Bound on i?Q 2 Vf iZjj 2 : The following elementary operator inequality is very useful 
to estimate the electron-electron interactions. 

Lemma 3.2. 

(41) ¥ -^<2(-A j -A k ). 

Proof. The unitary transformation induced by the following orthogonal transfor- 
mation of R 3 x M 3 , 

commutes with the Laplacian, and transforms \rj — r^\" 2 into 2 1 1 s | 2 . By Hardy's 
inequality, 

J- < -2A S < -2(A S + A t ), 
2\s\ 2 

and transforming back to the (rj, ^-coordinates yields (|41|l . QED 
We can then estimate: 



v i ) JV(JV-l) v i_ 



< TV (iV - 1) ^(-Ai - A,) = iV(iv - l) 2 ^(-Ai) 

< ^-1)^(4 + 1!^), 



and therefore, by similar arguments as before, 

VRo'V^y/Ro' < 2N(N-1) 2 ^/r q ~(T ± + ^(M + N))^r1 ) 
< N{N- l) 2 (M + iV + 2), 

on Ran (n_i_). 

Bound on VRq^ 2 : By the general identity ||AA*|| = ||^4|| 2 , we have: 

||V< 2 || 2 = WR^R^W < 2(||< 2 VX /2 H + HflS 7 2 V e 2 < 2 H) 
< (8Z 2 N + 2N(N-l) 2 )(M + N + 2). 
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Bound on VeS^R 1 ' 2 '- (Remember that we have shown that R > 0, so its square 
root is well-defined.) We first estimate ||ViZV||, as follows. Recalling the non- 
interacting resolvent i?Ni introduced above, we have that: 

< VRV < VRmV = VRl' 2 (l + R^VnR^ 2 ) * Rl' 2 V. 



Hence its norm can be estimated by: 



IV-RVII ■ 



\VR, 



1/2 || 2 

'0 ! 



1 - \\Rl /2 V n Rl ,2 \ 



from which we obtain an estimate for HVi? 1 ' 2 !! by taking square roots. Therefore, 
if B > S|2gJ = 46o as above, 

||V cff ,_Li? 1/2 || < ||Vi? 1/2 || < ^(16Z 2 N + 4N{N - l) 2 ) (M + N + 2). 



We now come to the proof of Q38JI. By Feschbach's formula, we have 



(H-0-'-RZ®R(0\ 



-R^V eS}± R 

— 1 1 -R c fF V e ff, _L 1 1 + ||-RV_L ) effi? ffV c ff ) _L-R|| 

(43)< Hi^ll \\V cS , ± R^\\ \\R^\\ + \\R™\\ \\V«,±&'*\\* \\R 1/2 \\ 2 , 

where we have used the following elementary estimate for the norm of matrices of 
operators: 

B )||<N| + ||5||. 
Hence if B > -BgSJ an< ^ ^ £ 4- <7 (-^rff)i we obtain 

\\( H BWffiB^ii < 1 ( ^ VRl/2 W i nvR 1/2 \\ 2 \ 

\\(H-0 -R cS ®R(0\\ < j= + B j 

< ,w/i n i/-> i c o 



with 4 = (32Z 2 N + 8N(N - 1) 2 )(M + N + 2). QED 

Corollary 3.3. (of the proof,) Theorem }!!. 1\ also holds, if we replace H and h e s+W 
by their fermionized versions Hf, /i e ff,f + VVf, and d^g(£) by the distance of £ to the 
spectrum of h e s r f + Wf . 

Proof. Simply write down the Feschbach's formula (|35|1 for Hi with respect to the 
decomposition / = n o fff +IIj_,£ of Hf, and estimate as in Ij43|l . where all operators 
will now have a sub-index 'f. Next use that P AS commutes with everything, and 
trivially estimate \\A { \\ = \\P AS A\\ < \\A\\, for A = R, R 1 / 2 and V ctf ^R^, except 
for \\R^g f ||, which will be estimated by 1 over the distance of £ to the spectrum of 
H cS , f + W f . QED 

The proof shows that in the fermionic case, Theorem 13 . II will at least be true with 
the same constants as for the boltzonic case. The optimal constants for fermions 
might be smaller, though. 

Remark 3.4. In the proof of theorerr f3.1l we systematically used Hardy's inequality. 
Alternatively, one can use, at least when N = 1, the bounds on the matrix elements 
of the Coulomb potential with respect to the Landau levels which were derived in 
[FW] , 
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4. Effective potentials for large fields 

The operator h cS + W = h^ M + W B ' M of Theorem IO acts on Ran(nf f j M ) = 
IQ M (W), aHilbert space which depends on both B and M, and which is canonically 
isomorphic to the space of i*^ -valued L 2 -functions on R N , 

(44) Ran(nf ff M ) = L 2 (K^). 

(Recall that F$ = Span {X B : m e S(M)}.) We will mostly suppress the Mi- 
dependence from our notations, M being fixed in our analysis. The potential term 

(45) v B s = nf ff vnf a , 

can be interpreted as an operator valued function of z = {z\, • • • , Zn) £ M. N , with 
values in the space of linear operators on F-S and acting in the natural way on 
L 2 (Mf ; F^).To get rid of the B-dependence of our Hilbert spaces we do a unitary 
re-scaling. Let us pose x = (xi, ■ ■ ■ ,xn) € R and similarly for y and z. Define a 
unitary operator U B on Ti by: 

(46) U*rf>(x, y, z) = B N / 2 i>(VBx, ^By, z). 
Since X B (x, y) = B N ^X^(VBx, y/By), it follows that 

TrB*yrB ttB _ ttI 

u xy L1 eS u xy ~ il eff- 

Let us write V B S in multi-particle form: 

(47) V C B S = -J2zvf+J2vf k , 

3 j<k 

with Vf and Vf k defined by GO}. Then U%VfU* = I^ B U^\r J \- 1 U^ s = 
y/BVl{y/Bzj), with 

(48) v}{z) = i^ ff 1 =n^ ff , 

/x 2 + y 2 + z 2 



and likewise for TAf : U^VgU^z) = y/BV^VB^Zj - **))> with 

1 

^/(zj - a; fc ) 2 + (?/,■ - y fc ) 2 -I- z 2 



(49) y& (z) = n^, 1 =n^ ff . 



The operator 

(50) /i efI := U x *H eS U B y , 

will now act on the fixed, B-independent, Hilbert space, L 2 (R N ,F^), and 
h B s = -~A* - ^ zVBV^VBzj) + VbV^(Vb( Zj - z k )) 

3 3<k 

= -±A z + VBVt s (VBz). 

The next step will be to examine the asymptotic behavior of y/BV^ s (y/Bz) as 
B — ► oo. The main idea is contained in lemma |4~T1 below. We introduce the free 
Laplacian on M. N , 

(51) /too = -^A„ 
and its resolvent: 

(52) i?oo(-« 2 ) = (hoo+a 2 )- 1 . 
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We will need this resolvent both in dimension N and dimension 1. To distinguish 
between these two cases we will, in the 1-dimensional case, systematically use [3 2 
as spectral parameter instead of a 2 , reserving the latter for the multidimensional 
case. 

If u is a function or tempered distribution on WL N , with values in some auxiliary 
Hilbert space F, then || -Roo ( — a2 ) s u \\l 2 (m n ;F) is a norm on the s-th Sobolev space 
H S (R N ;F). A linear operator A sends H S (R N ;F) continuously into H~ S (M. N ; F) iff 
the L 2 -operator norm \\Roo(—a 2 ) s ^ 2 ARoo(— a 2 ) s ^ 2 \\ is finite. The case of interest 
for us will be s — 1. We will also need the Fourier transform J 7 , but only in 
dimension 1, for which we normalize as follows: 

■F(u)(C) = / u{z)e- tz ^dz. 
Js. 

There will consequently be a factor of (2tt)~ 1 in the inversion formula. 
Recall that 

Pf(^) =^(sgn(x)log|x|), 

with the derivative in distribution sense. Let F be a finite dimensional complex 
Hilbert space, and L(F) the space of linear operators on F. 

Lemma 4.1. Let be an L(F)-valued tempered distribution on R, such that its 
Fourier transform can be identified with a locally integrable function TX) — J-v(£). 
Assume also: 

(i) There exist Co, C\ € L(F) and a > 1/2, such that: 

(53) ^o(0 = -Colog|C|+Ci + 0(|C| ), C^O. 



(H) If 

e(C) :=^t>(C)+Cblog|C|-Ci, 
denotes the error in the approximation ^53j) . then 

(54) 6 := / d( < oo. 



For each A > let 



(55) Ooo.a :=C logA-J+ic -Pf ^ J + ( 7 C + d) • S, 

where S is Dirac's delta- distribution in 0, and 7 = T'(l) is the Euler constant. If 
Roo{—/3 2 ) denotes the free resolvent in dimension 1 and [3 > 0, then 

(56) ||i? 00 (-/3 2 ) 1/2 (Ao(A-) - 0oo , A ) i? 00 (-/3 2 ) 1 /2|| < t!^. 

Remark. Observe that the integral l|54|) converges in 0, since we assumed that 
a > 1/2 in (ESJ. 

Proof. It is known that 

(57) ^(loglClH-Wrjr) -7*0, 
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where 7 is Euler's constant: cf. e.g. Schw ; . Therefore 

tw = Co(logA) S-CoT- l Oog\C\) + Ci6 
= ^- 1 (-C log(|C|/A) + C 1 ), 

and 

T (i?oo(-/? 2 ) 1/2 (Ao(A-) - tw A ) R 00 (-f3 2 ) 1/2 ) T~ x 
is an integral operator with kernel: 

( 58 ) ^Z ^2/o I /32M/2 e 



2^(c 2 /2 + /? 2 ) 1 / 2 v a ; (c'Va + z? 2 ) 1 / 2 ' 

since multiplication by a distribution a(x) becomes an integral operator with kernel 
(2Tv)~ 1 J 7 a(<^ — (') after conjugation by T . Since conjugation by the Fourier trans- 
form does not change the operator norm, it follows that the norm in (|56|l can be 
bounded by the Hilbert-Schmidt norm of (|58|l . whose square equals: 



TV 



V2 fMv/W^^f Hf C1V2 



7T/3 J R T? 2 + 8/3 2 ' 7T/3A y R r? 2 + (/? 2 /8A 2 ) ' " tt/3A ' 
Here we have used the elementary integral identity: 

where a, b > 0. This finishes the proof of lemma fOl QED. 

We will apply the previous lemma to our potentials fTT)| , but before doing so we 
first state and prove a weaker variant, which will be used to prove Theorem 11.11 
Let us introduce the (numerical) constant: 

f 1 f (l lQ gH l + 2 ) 2 * V /2 

Numerical evaluation of the integral (using either Mathematica or Maple 8) gives 
Lemma 4.2. Let t> = t>(z) be an L(F)-valued tempered distribution on M such that 
(61) C := sup || (| log (CI I + I)" 1 J*(C) II < °°- 



Then for all A > e and all e > 0, 

A 

.log A 



(62) || (sA + e- 1 )- 112 i-tKA^) (-eA + e" 1 )- 1 / 2 || < C m C v (| loge| + 2) . 



Proof. Conjugating as before by the Fourier transform, and estimating the operator 
norm by the Hilbert-Schmidt one, we find that the left hand side of (|62|l is bounded, 
by the square root of 

1 .„\ 1 „ .,1s 1,2 



^ 1571 can also easily be shown directly, using the observation that Pf and — 2jF _1 (log |£|) 

are both solutions of the distributional equation xA = — sgn x and therefore only differ by a 
multiple of 5, which can then be computed to be —2-y. 
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By (|6*l"1) we can bound 



1 "^/A)ll 2 < c?fi^!?li + i + . 1 



(logA) 2 " v " y " " D V lo S A lo S^ 
< C D 2 (|logH | + 2) 2 , 

since we suppose that A > e. Hence, using Ij59(l again, we find that our norm is 
bounded by the square root of 



Cl f (| logM | + 2) 2 Cl f (| logM | + |log £ |+2) 2 

C 2 ,2/ *? /■ (| logM I + 2) 2 , _ ^ /(loge) 2 



by (jnOJ). Since C|^J > 1/2, we see that will be bounded by CjgrjC (| loge| + 2), 
as claimed. QED 

The next step will be to apply lemma 14.11 to the potentials Vj~ and Vk , with 
A = \pB. We introduce the independent tempered distribution q = q B , and linear 
operators Cf, C] k € L(F£) by: 

(63) g B (z)=logi?^(z) + Pf ^ 

(64) Cy := C- M := -^ ff log (I (x) + y 2 )) U] s , 



(65) q fc := C*f := -H^log (!((*,• - x k f + (y 3 - Vk f)) U\ s ; 

Observe that (|6~4l and llB5|) are related to |Q and l(T3|) by conjugation by U^ y . See 
also remark f 1.41 for a physical interpretation of these three terms. 

Lemma 4.3. Let R 00 (—f3 2 ) — (—^A z + /3 2 ) -1 , (3 > be the free resolvent in di- 
mension 1. There exists a positive constant Cjg^J := C fi^ (M) > only depending 
on M, such that for all B,(3 > 0, 
(66) 



||i?oo(-/? 2 ) 1/2 (VBV^VBz) - (q B (z) + CJ S(z))) i?oo(-/3 2 ) 1/2 || < 

and 
(67) 



j j kjuv m ; 11 ^ 7^ 

i/ie norm being the operator norm on L 2 (R, i^). 



(-z? 2 ) 1 / 2 (^i(^) - ( g B (z) + *(*))) i? fl0 (-/? 2 ) 1/2 || < 



To simplify future estimates, we have taken the same constant in both inequalities. 



Proof. Recall the formulas JIB} and (gHJ for VHz) and V^z). We need the 
asymptotics of their Fourier-transforms at 0. By |AS1 9.6.21], the Fourier transform 
of (1 + z 2 ) -1 / 2 equals 

^((l + z 2 )- 1 / 2 ) (C) = 2^o(|C|), 
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where Kq is the Macdonald function. Since the projector H* s effectively only acts 
in the x and y- variables, it follows that 



= 2ni ff ^o( v /^ + y 3 2 -C)n c 1 ff , 

with a similar formula for J~Vp.. 
Now it is known that 

7,o(|CI) = -log|CI+log2- 7 + 0(||C 2 log|C| |), ICI-0, 

and that -KodCI ) is bounded on |£| > 1 (even exponentially decreasing there): see 
e.g. |ASI 9.6.13]. It then easily follows that, as ( — ► and as operators on Ran nj ff , 

(68) ^ 1 (C)^-21og|C|-2 7 + q i , 

^(C)^-2io g |c|-2 7 + q fe , 

with an error of 0(|£ 2 log|£| |). An appeal to lemma ITTl with A = y/~B and with 
Co = 2 and C\ = — 2 7 + C™ respectively C\ = — 2 7 + Cj k , then finishes the proof. 
QED 

Remark 4.4. An explicit computation of the matrices of C™ and C? fe with respect 
to the natural basis X^m 6 S(M) shows that C" and Cj k do depend on their 
indices j and j, k, respectively. 

We will likewise need lemma for = Vj . We can without loss of generality 
assume that j = 1, by permutational symmetry of S(M). As we have seen above, 
TVl{C) = 2H 1 cS K {\CW4+yl)^ls- It can easily be verified that || (| log |C| I + 
l) _1 .Ko(C) lloo = 1) so that, for example, 



(69) C v i < C m := 2 + 2|| 



log J x\ + y\ 



The operator norm on the right and side can be evaluated explicitly, and behaves 
asymptotically for large positive M as 21og(M). 

We next extend lemma |4~51 to multi-particle potentials. Let us define the multi- 
particle potential vc by 

(70) v c (z) = vS(z) = -Zj2(q B (z j ) + Cfd(z j )) 

j 

j<k 

Lemma 4.5. Let Roo(~a 2 ) = (— |A 2 + a 2 ) -1 , a > 0, be the resolvent of the free 
Hamiltonian in M. N . Then 

1 1 - "uu v — / I * ~ ' r en \ ' / -ks \- j t ~ "uu v — / ! I - r 

v ' v< 

where 

(72) := Cffl (N, Z, M) := Cfffl N 1 '* \Z + \{N - 1) 

Proof. We split both potentials into their 'electron-nucleus' and 'electron-electron' 
parts: 

V 1 ^ = Vff + V 1 ^ 



(71) || J R o(-a 2 ) 1/2 (VBVyVBz) - « c (z)) i? 00 (-a 2 ) 1/2 || < 
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and similarly for v C - v c = v c , n +vc,e = v c,n+ v c,e- Writing V B v {z) for \fBVl v {\/Bz) 
(with a mild abuse of notation), where v — n or e, we bound the left hand side of 
J7T} by 
(73) 

||i?oo(-« 2 ) 1/2 (V^n-vcn) R m (-a 2 ) 1/2 \\ + \\Roo(-a 2 )^ 2 {V c % e -v c ,e) Roo(-a 2 n 

and estimate the two terms separately. Let Roqj(—(3 2 ) be the 1-dimensional resol- 
vent in the variable Zj, with a (3 which will be picked below. We will simply write 
i?oo for Roo(— & 2 ) and Roo,j for Roo,j(—f3 2 )- If we put 

AV j :=V 1 B {z j )-q B {z j )-Cn{z j ), 



and 



then, by ^EJ, 



C M Z y nl/2n-l pl/2 _ C M Z p / ^2/ ! A , ^ 



< 



3 

CfiftZ |C| 2 /2 + iV/3 2 _ C m Zm^(l,N(3 2 /a 2 ) _ C m ZN^ 

v/TJS 1 /*^ |£| 2 /2 + a 2 y%B 1/4 V^^ 1/4 

if we pick (3 = a/yN; this choice actually minimizes /3 -1 / 2 max(l, N(3 2 /a 2 ) as a 

1 /2 1/2 

function of /3 > 0, as is easily checked. Similar estimates show that — R 0Q AV n R QQ 
is bounded from above, in operator sense, by the same number, and we therefore 
conclude that the first norm in J73J) is bounded by CjftftZ W 1 - /4 / 'y/aB 1 • 74 . 

To estimate the second term of (1731) . we will use the following lemma, which is 
analogous to lemma l3~2l from section 2. Let Aj = —d 2 /dz 2 , A k = —d 2 /dz 2 . 

Lemma 4.6. Let t> = t>(z) be an L(F)-valued distribution on R (F a finite- 
dimensional Hilbert space), such that Rqo{— /3 2 ) 1//2 t)i£oo( — ft 2 ) 1 ! 2 is self-adjoint, for 
(3 > 0. Then for all [i > 0, 

(74) || (—A, - \A k + m 2 )- 1 / 2 *(z 3 - z k ) {~ l -A 3 - ^A fe + /i 2 )" 1 / 2 || 

< I||(-I A . + £)-V» ,(.)(- 1^ + ^-1/2 11, 

where the norm on the left hand side is of course taken in L 2 (M. 2 ,F). 

Proof. We use a similar change of variables as in the proof of lemma 13.21 s = 
(zj — 2fc)/V2, t = (zj + Zk)/V2. Then, with ~ denoting unitary equivalence, 

1 1 \~ 1/2 / 1 1 ^~ 1/2 

' ' ' • 2 \ „/ „ s / A 1 A i ,,2 



--A J -~A fe + /^ J n[zj-z k ) \—-Aj — -Afe + fi' 



11 \ ~^ 2 / 1 1 x ~^ 2 

-2A s --A t + M 2 J o(V2s) (--A s --A t +, 

i/i i 2\ - V2 / -, -, 2\ - 1 / 2 

1 ( 1 a 1 a , M 2 \ .., x Ma 1 a , M ^ 



2V"2 A ^2 At + Y ° (S) -2 A '^2 At 



Observing that 



1 u 2 \ 1/2 / 1 1 M 2 V 1/2 

2 2/ \ 2 2 2 1 
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on L 2 (R 2 ), the lemma follows. QED 
Let us write 



(T.-») /?„„.,■/, : /?o..., ( 7.(- /ri = ( -i(A, +A fc )+ At 2 



-1/2 



the 2-dimensional free resolvent, where will be optimized at the end of the proof. 
Recall that Rqq = R^i—a 2 ), and put 

AV ]k = Vf k ( Zj - z k ) - q B ( Zj - z k ) - q k S( Zj - z k ). 

Then, using lemmas IQ1 and ET31 R 1 ^ 2 (V B fie — v^^R^q 2 can be estimated from 
above as follows: 



^OO 2 X] AV 3 k ^OO 2 ~ X! (^OO^OOJfe) ( R 00 2 jk AV jk R 00Jk) ( R 00,jl R 

j<k 

^El V^nl/2p-l pl/2 

23/4^^1/4 ^00 ^OOjfc^OO 



/2 
00 



j<k j<k 

a 



2\ 



= 23/4^/4 2 A - + ^-j 

^ r.-J//l nl //I ' ~i=- maX i ' T7~^" S ~ 



where we minimized the right hand side over fi by choosing [i = a^/2/N . The 
similar upper bound for-i? 1 / 2 (V e - v^^R 1 ^ 2 gives the desired estimate for the 
second norm in l|73|) . and combining the two estimates, we have proved lemma 14*31 
QED 

We now derive a similar estimate for W B as B — > oo. 

Lemma 4.7. Lei i?oo( — c* 2 ) & e /ree resolvent in dimension N , and let U = U B 
be the unitary transformation defined by \4b\ - Then, if £ < 0, 

(76) ||i?oo(-« 2 ) 1/2 ^W s f/i?oo(-« 2 ) 1/2 || < 

a\J B 

with 

(77) C^:=C^ (N, Z) := 2^ 2 N^ 2 (z 2 + 

Proof. Recall that W B = -V eS ,±RVx, e &, where i? = i?(£) = (H ± - on 
Ran Uf. Hence 

J7*W B J7 = -BV^ ± (VBz)U*RUVi ieS (y/Bz), 

where ± = n^VH^, and similarly for V]_ cff . Since for £ < 0, < R < 2/B 
(see section 3), we have that, letting V(-, z) be the function (a;, y) — > V(x, y, z), 

o<-u*w B u < 2nJ ff v(-, Vbz)u±v(-, Vsz)n^ ff 
< 2n 1 cH v(-,VBz) 2 n 1 cfl 

(78) < 4(nlM-,VB Z ) 2 nl s + nl s v e (-,VBz) 2 nl s y 
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As in the previous lemma, we treat the two terms separately. By Cauchy-Schwarz, 



V n (x,y,^z) 2 < Z 2 NJ2^- 



where p 2 = x 2 +y 2 . As before, let R 00 .j(—P 2 ) be the resolvent of h 00 .j = -(l/2)d 2 /dz 2 
on R. We first estimate the L 2 -norm of each 

(79) R oA-P 2 ) 1/2 nl s (p 2 3 + S^ 2 )" 1 n o 1 ff i?oo J (-/3 2 ) 1/2 , 

by conjugating with the Fourier transform T . Since T ((1 + z 2 ) -1 ) (C) = 7re _ ^l, 
(|79|l then becomes an integral operator with kernel 

(80) \b~^ 2 (f + P 2 )"' 2 (n^ pfe-^-^-nA n ^ (± + /3 2 )" /2 . 

The norm of H80JI can be estimated by its Hilbert-Schmidt norm, whose square can 
be bounded by: 

ilin^-n^H 2 ( I (£ + / 3 2 r 1 ^ 2 - 



4B" cn p 3 on " V7r 7 2/3 2 B ' 

where we have put 

(81) c m := || Il^nk ||. 

Note that C jg^j is independent of j, because of the permutational symmetry of 
E(M). It follows that: 

(82) i?oo,,(-/? 2 ) 1/2 (p 2 + Bz^Ul, RooA~P 2 ) 1/2 < 
and therefore 

< i? 00 (-a 2 ) 1/2 ni ff V„(-,VB) 2 ni ff Rooi-a 2 ) 1 / 2 

iS3% J&BV* R ° 0i ~ a )E(V J +/? 2 ) = ^ 1/2 ^oo(-^ 2 )(feoo + ^ 2 ) 
< ,C m Z 2 N^ 2 



V^aB 1 / 2 



if we choose /3 = a/yN. The same inequality then holds for the norm, since the 
operator we estimate is positive. 

We next treat the interaction term i?oo(tt 2 ) 1//2 V e (-, V^z^-Roo^ 2 ) 1 / 2 in a similar 
way as in the proof of lemma Ef31 First, by Cauchy-Schwarz again, 

where we have put p 2 k = (xj — Xk) 2 + {yj — Dk) 2 ■ Since the rotation (rj, rk) — > 
(2 _1 / 2 (r :) — rfe), 2 _1//2 (r 3 +7"^)) commutes with IIg ff (since it commutes with Hq 
and with L z ), we find, after a unitary transformation, that 

II i?oo Jfc (-^ 2 ) 1/2 n^ (p 2 k + B( Zj - zkfy 1 n c 1 ffJ R o Jfe (-/i 2 ) 1/2 II 
< II UooA-^) 1/2 K s (p 2 + Bz 2 y x nyw--^ /2 11 



< 



2/^5 ' 
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compare the proof of lemma WM Hence, using similar arguments as before, 
i?oo(-a 2 ) 1/2 Kg V e (-,^) 2 Ill s i?oo(-« 2 ^/ 2 



*C m N(N - 1) 

< — , 7= Roo(-a ) > (/loojfe + A 4 ) 

4/iV-B ^ 

.C^A^-l) 2 JV/A 
- Wl )(*» + —) 



< 



if we choose /j, = a^j2/N . Adding this estimate to the one for V 2 , and remem- 
bering the factor 4 from l|78|) , we have proved l|76l ) with the constant C f^ = 

2 3 / 2 7rCj^jA r3 / 2 (z 2 + (jv ~ 1)2 ) . A priori, CjgJJ might still depend on M, but it in 
fact does not, as we will finally show. We compute C |g*[l m the Landau basis l|28|) 
of Fm 1 , with respect to which n^ ff p^ 1 LT^ diagonalizes: 

ClRTl = max ^ r / pi m e- p * /2 d Pl 

J l o<m<M2 m m! /« 1 

1 



3G 



0<m<M mlv^ 7o 

r(m + |) 



o<m<M^2r(m + l)' 



1^ _ 1-3-5 (2m- 1) /l 



It is known that 



2) 2 m \2 y 

cf. e.g. |AS| . formula (6.1.12), page 255. Using this, one easily finds that 



which completes the proof of the lemma. QED 

To prove Theorem ll.5l we will need to control the Sobolev norm of V® — vf . This 
is done in the following lemma, which we formulate in slightly greater generality 
than needed, with an eye to future applications. 

Lemma 4.8. Let V^ 3 := Vj 3 ' be defined as in \1U(I . c > and let a c — a c (B) be 
the unique positive solution of 



(84) ac = ri og i^ , 

c V a c J 

Then the constant 
(85) 

crm ■■= K + 9io g 2 (2) + ^ + sup 1^(0 + 2io g (ici)| 2 + ^ sup ijn/^OP 

*** \ 71 \C\<i n Kl>i 

is finite and depends only on M. Moreover for all c > and all B > 

(86) ||i?oo(-« 2 )^(n B -ca c (S)^i?oo(-« 2 )^|| < 



1/2 
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Observe that if c = 2, then a c (B) is the a(B) from Theorem ll.il cf. I©- Also note 
that the constant C |gg| is independent of c > 0. 

Proof. The first statement follows at once using the explicit knowledge of TV^ 
obtained in the proof of lemma 14.31 To prove the estimate (|86[1 we introduce the 
auxiliary function X by 

FVfiO-ca^iO = ^ 1 1 (^) + 21og^-21og^ 

VB VB u c 



VB a c 



where in the first line we used equation l|84Jl . The Fourier transform of Vj 3 — 
ca c (B)S acts as convolution by a function, and we estimate the norm of Y := 
Roo{— a c) 3 (^i B — ca cS)Roo(— a?) 5 by its Hilbert-Schmidt norm, as in the proofs of 
lemmas 1001 It follows that 

||Ff <V2 / |^(C)-ca c ^(C)l 2 dC < 2V2 r |A(^)| 2 + |21ogg| 2 ^ 



First, 



8«2 4 


-c 2 




L 8o% + C 2 


log 


jC[ 







8« 2 + C 



1 '• ' 



, s dC=- 7 =- T (^ 2 + 91og 2 (2)). 

We next look at the contribution of X , which we split in two parts: 

00 l*(jg)l 2 <*C j_ r°° |x(Q| 2 dC J_ /'°° |^(C)l 2 +4(iog|CI) 2 , f 

8a2 + C 2 85-1^ + C 2 " A C 2 

< 2 ( sup|jry(c) |2 
vB V I CI > i 




since 



and 



|£|2 aC » Z ' 



8a c + C 2 |C|<i Jo 8o^ + C 2 

< * sup|X(C)| 2 . 
4y2a c |^|<i 

The rest is now elementary. Notice in particular that sup B>0 a c (B) / \/~B = 1. QED 
The limit potential fTTHl suggests defining an effective Hamiltonian he — h^ by: 
(87) he '•= h 00 + v c - 

As it stands, this is just a formal expression, and our first task is to give a meaning 
to he as a self-adjoint operator on L 2 (R N ; F^) . We do this by showing that vc 
is form-bounded with respect hoo, with zero relative form-bound. Let (•, •) denote 
the duality between distributions and test functions. 

Lemma 4.9. The quadratic forms u — > (S, \u\ 2 ) = |w(0)| 2 and u — > (Pf(l/|z|), |m| 2 ) 
are well-defined on H (K), and form-bounded with respect to hoo, with relative 
bound zero. More precisely, we have for all e > that 

S< |(-eA, + E- 1 ), 
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and 

(89) PfCN" 1 ) < C m (\ loge\ + l)(-eA 2 +6" 1 ) 



wh 



ere 



(90) :=A/^ + 2(log2)2 +7 . 



r 2 
2 

Proof. This is well-known for 8. For Pf ( | ^ | 1 ) we first note that, since Pf (|£| _1 ) = 
— 2T~ 1 (\og ICI) — 27^0, it suffices to prove the form-boundedness of .F _1 (log |£|). 
The latter will follow from: 
(91) 

|| (sA + s^)-^-\\og \C\)(-eA + e-')-VY < \^gsf + + 

for all £ > 0. To prove (|91|l . observe that after conjugation by the Fourier transform 
J 7 , and estimating the operator norm by the Hilbert-Schmidt norm, the square of 
(191(1 can be bounded by: 

1 r r (losic — /?i) 2 

(92) (^LL(ee + e-i)(e V i + e-i) d ^ 
Changing variables and using l|59(l . we find that i|92|) can be bounded by 

J_ f (log ICI) 2 , c = J_ f (logical) 2 ,, 
2 7 r7 R eC 2 +4e- 1 ^ 2n J R C 2 + 4 ? 

dC If (logC) 2 _ l n ,2 , 7T 2 +4(l0g2) 2 



I" A C + 4 TT J R C +4 2 

Hence, using iJSSJ, 

|(pf(N- 1 ),| M | 2 )| = 2. ((^(logic-D.K 2 



8 



1/2 

2 




1 ^ ( 2(loge)^ + y + 2(log2) 

which implies ijEHJ since v^/Cjg^J < 1. QED 

Because of the terms C™<5(zj) and Cj k 8(zj — 2^.) we have to extend the first part 
of lemma R~9l to the vector- valued case, but this is immediate: we just note that if 
C a linear operator on F^ (or any finite-dimensional vector space, for that matter), 
then the right interpretation of C8 as quadratic form on ff 1 (K,F I J [ ) is given by 
(8(z), (Cu(z),u(z))), where (•,•) is the inner product on F^. 

Finally, we lift lemma fOl to R w . Recall that if L : M. N — * R is a linear map, then 
the pull-back L* A of a distribution A on R is well-defined, and can be computed by 
going to linear coordinates z' with respect to which L(z') = z[. It then immediately 
follows from lemma FOl that L*8 and L*q B will be form-bounded with respect to 
Hqq on 1*, with relative bound 0. Taking L(z) = Lj(z) = Zj and L(z) = Ljk(z) = 
Zj — Zk, we see that the sesqui- linear form t B {u) given by 

(93) tg(u) = ±\\Vu\\ 2 -Zj2(L*q B , \u\ 2 ) + ((L*8, (C>,u)) 

j 



j<k 



(L* k q B , \u\ 2 ) + (L* k 8, (CJ k u,u)) 



is well-defined on H 1 (R N ,F^), and bounded from below by — C||m|| 2 , for some 
constant C, depending on B, Z, N and M. By the Kato-Lax-Lions-Milgram-Nelson 
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Theorem (cf. e.g. jEHI Theorem X.17]), t% defines a unique self-adjoint operator, 
which we will call he — h^' U , and informally write as (|87J) . In Appendix A we will 
give a characterization of the operator domain of . Similar arguments will define 
hf as a self-adjoint operator, cf. |BD| . 

5. Proof of Theorem ll.il 

We will first compare the resolvents of h c g :— h^ M and of /i e ff + W, where 
W := W S,M (£) was defined at the beginning of section 3. Put 

r cS := r cS (0 := (h cS - _1 , B™ ■= B™(0 ~ (h cS + W - , 



AC 



'fen 



2 ' 

a[Cr 



and let 

(94) B M 

where a — a(B) is the function defined by © , and 

(95) ^ : = 2 {^ + 1 ) C m 
with £ e ff = e e s(Z, M) the unique positive solution of 

(96) ZC m C vl e(\loge\ + 2) = ± 

where C v i is defined by (|6I[) with t> = V± (see l|69[) for an upper bound). Note that 
both constants only depend on N, Z and M, and this in a controlled way. 
We then have: 

Theorem 5.1. If B > Sj^J , £ < and 

(97) c ce ^= < d cB (0 < ia 2 , 



with a = a(B) is as in Theorem then £ G p(h c s + W), and ||-Rrff(£)| < 

2||r c ff (£)|| . Furthermore, 

(98) ll^(0-reff(0ll<ceff- " 



Proof of Theorem \5.1\ It clearly suffices to establish Theorem 15. II after conjugation 
by U^y, defined by |@EJ- To simplify notations, we will simply denote the conjugated 
operators by the same letters as the original ones. Using the symmetrized resolvent 
formula we estimate 

Wf« -_^MI^ 1 ( ll*eff(OII 



where 

(100) *T eff (0 := IreffCOI^WCOr-eff^) 172 , 

with the convention that A 1 / 2 := sgn (^l)IAI 1 / 2 , if A is a self-adjoint operator 4 . 
Now let [L < inf (cr(/i c ff), to be specified later. The following elementary lemma will 
allow us to replace £ by /i in (|100|) . 



note that since £ is not necessarily below the infimum of the spectrum, r e g(£) is not necessarily 
positive, and we use the symmetrized resolvent formula in the following form: R^L = r^^(l + 
Ireffl^Wr^VVeffl 1 ^ 
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Lemma 5.2. If p < inf <r( ft» e ff)> then for all real ^ in the resolvent set p(h c g), 
(101) || r eff (0 (&eff-A0 II <max(-^J-,l 



Proof. We distinguish two cases: inf c(/i e g) < £ < and £ < inf cr(/i e ff)- (Observe 
that inf a(/i c ff) < 0, by the HVZ Theorem, since this is already the case for N = 1). 
In the first case, let (£_ , £+) be the largest open interval in p(h c g) which contains £. 
Since [0, oo) is in the spectrum (it is already in the essential spectrum), £ + < 0. It is 
easy to see that the function x — > | a; — yuj / 1 x — £ | is increasing on ( — oo , ) n <r ( /i e ff ) 
and decreasing on (£ + , oo) n o~(h e g). It follows that 

||r eff (0(^-/,)||= sup !^=rnaxfl^i ! l^i > )< ^ 



as was to be shown. One shows in a similarly way that (|1(J1[I equals (inf a(h e s) — 
/Lt) / (inf a(h e g) — £) < |^|/d c ff (£), if /i < £; < inf cr(/i ff ), and is equal to 1 if £ < //. 
QED 

Substituting Id = r ef r (^i) 1 / 2 (/i of f — /i) 1 / 2 = (/i c ff — m) 1 ^ 2 ^^) 1 ^ 2 at the appro- 
priate places in formula (|100|l . we see that if £ < 0, 

(102) \\K cS (0\\<^(-r4^A) \\K eS (»;Ol 



d eS (0 ' 
with 

K<b(ji;S) := r cff (/i) 1/2 W(Or c ff(/i) 1/2 - 

Repeating the same argument for K e s(p;^) using Id = ((/loo + a 2 ) 1 ^ 2 -Roo(^a 2 ) 1 / 2 , 
we obtain from lemma l4~7l that 

(103) \\K(n;0\\ < ^||r cff (/i) 1 / 2 (/ l0 o + « 2 ) 1/2 || 2 . 

We will now estimate the norm on the right hand side, for suitably chosen p. 
Lemma 5.3. Assume B > e 2 . Define 

N 

cff 

where a is as in Theorem M.lY and where e — e c s is the unique positive solution to 
the equation 1.961) . Then p c g < inf (j(ft. ff ), and 

(105) ||r cf f(Mcff) 1/2 (/ioo + « 2 ) 1/2 || < V2. 

Assuming the lemma for the moment, we continue with the proof of Theorem 
O we have, by JUSJ), JTHJ) and (fT77s|) . that if d e & (£) > c # a/VB~, then 



(104) f, eS =p efi (N,Z,M) :=-— I + 1 



K .„ (0II < Jm „^ ( ^ + 1 )_J i _ il }^ fi i, 



provided that c# satisfies 



c *- 2 (4 + 1 ) c El'» d 



Since a{B)\/~B is increasing, and since a(B)\^B = x iff £> = x 2 /(4a(a;/4) 2 ), the 
last inequality is implied by B > BjT^J- Choosing c# = c G ff defined by Ij95(l . we 
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conclude that if £ is such that d e s(0 > c c ga/VB, then by ij^H) , i|102l) . I|103|) and 
our choice of \i = /i e s, 



R%-r eS (0\\ < 



a 2 ( N 2 \ 1 1 a 



4ff(£) I 2 V2e e s / 4ff(0' J aVB 



: 2(7^,(^ + 1^ 



2£ 2 ff / dcff(£)W 

provided that d c ff (£) < a 2 /2. This proves Theorem 15. II modulo that of lemma l5~31 
QED 

Proof of lemma \5.'A We will use a scaling argument. If we let ~ denote the 
unitary equivalence induced by the change of variables z — > z/a, where a > is 
for the moment a free parameter, and if we write Vjj for ^/BV^, s (VB-), then if 
fi < inf o-(h cS ), 

< (^oo + « 2 ) 1/2 rcff(M)(^oo + a 2 ) 1/2 

- (V, + 1) 1/2 (-^A, + «' 2 Vf ff (-) - o- 2 M ) ~* (^oo + 1) 1/2 

(106) < (hoo + iy/^-^+a-^J-)-^ 2 ^ 1 (ftoo + l) 1 / 2 , 

Vg ff n being the attractive part of V c ff ■ We now choose a = log(y/B/a), as in @, 
and put A := \[Bja. Notice that B > e 2 implies that A > e. Then, by lemma FOl 

(VB \ „A A 



> -2C m ZC v i s(\ loge\+ 2) (-Ia z + ^\, s>0 



1, N 



2 4e e ff 

uniformly in A, if we choose e := e c ff = e c ff(Z. M) such that holds. We now 
take 

a 2 / N_ 

2 

cir 



Then it follows that /i e g < inf o~(h e g), since /i e ff — n s > ^(^oo + 1) > \- Further- 
more, (|106|) can be estimated by 

(^oo + 1) 1/2 (iftoo - (J- + ) ) 1 (*oo + 1) 1/2 = 2, 

which implies H105(l . QED 

Proof of Theorem It now suffices to combine Theorem 13.11 and Theorem 15.11 
while carefully keeping track of the constants. First of all, 197(1 . implies that 
ll-R^II < 2||r cff (£)||, and therefore d^(£) > d off (£)/2. Furthermore, £ < if 
d e s(0 > 0, and if 

(107) B>B cS :=ma,x{B m ,B^,e 2 }, 
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then 



< \\ {H - 0- 1 - R%(0 © (H ± - 0" 1 1| + - r eS (0|| 

9 

C H7\ + — ) , 7= ( usin g d eff(C) < a 2 / 2 a g ain ) 

2 



a(B cS )J d^YVB' 

since a(-B) -1 is a decreasing function of B. This proves Theorem l 1 . 1 1 wit h a constant 
Coff which is equal to 

(1^ ■= rrrm 4- 

a{B c g) 



6. Proof of Theorem IQl 

As a first step, we will compare the resolvents r e ff (£) of /i e ff and J~c(£) := (^c ~ 
0" 1 of ft c := /ig. Recall, that d c (0 := dist (£,o-(ft c ))- 

Theorem 6.1. Let a = a{B) be defined by 0). There exist (computable) constants 
B' c , C' c > 0, onZi/ depending on Z, N and M, swc/i t/iat /or aZZ £? > and o/Z rea/ 
£ < satisfying 

(108) dc(0 > Cj, a 3/2 £T 1/4 , 

</ien £ G p(H c ff), with ||r ff(£)|| < 2||rc(0ll- ^ n addition, letting 5 



C'i = maxiC'c, 4Cgg + l) } > C c 

then if 

(109) efc(0 > C£ o?l 2 B- x l\ 
we also have that \\rc(0\\ ^ 2||r e g(£)|| ■ 
Finally, if 

(110) C' c apPB-V* < d c (0 < \a\ 
then 

a 3/2 



(HI) \MZ)-r c (&)\\<C c 



dc(0 2 B 1/4 



Proof. The proof is similar to the proof of Theorem l5.ll with however some technical 
changes, due to the fact that vc is not homogeneous of degree —1, and that its 
electron-electron part is not positive anymore. As before, we conjugate all operators 
by U^y, keeping the same letters for the conjugated operators. 
Arguing as in the proof of Theorem 15. II one shows that 

ii <e\ <t\\\<r 1 \\ K c(i)\\ 
Fcff(f) - rciOW < 



^with e e ff defined in 1961 
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where Kc ■= |?*c(£)l (^eff — ^c) r c(0 1 ^ 2 - Using lemma |4"31 we find that, for any 
p < inf (r{h c ), 

(112) \\Kc{i)\\ < ^3}3= max|^L,l| \\r c (^ (h 00 + a 2 )^ ■ 

We then use the following analogue of lemma 1^1 of which we only state a qualita- 
tive version. 

Lemma 6.2. There exists a constant vc = vq{Z, N, M) > 1/2 such that if B > e. 
and if He : = -"c« 2 , ol defined by then pc < \aia{hc), and 

(113) ||rc(^c) 1/2 (^oo+a 2 ) 1 / 2 || 2 <2. 

Proof. As before, we will use scaling. However, contrary to 5(z), the distribution 
Pf(l/|z|) is not homogeneous of degree —1 on R. In fact, if p a is the dilation 
p a {z) = a~ 1 z on R, a > arbitrary, then the pullback of Pf (1/| • | by p a equals 

(114) piPf (ri) =aPf (n) - 2aloga5 - 

Let us split the potential vc '■= Vq of he as 

v c = log 5 v s + vq, 

where vs is defined in l|17|l and 



the (pseudo-) Coulombic part. If ~ denotes unitary equivalence with respect to the 
dilation p a (on M. N ), then in view of our choice of a 

2 / logS-21oga 1 _ 2 \ _ 2 ( _ ( 1^ 



h c ~ p a /loo H W5 H - a A* = a (^oo + 2w«5 H Wo 

\ a a J a 

> a 2 f ~/ioo - 6 - a^ 2 M 

form some 6 > depending only on Z, N and M, since by lemma 14.91 we know 
that 2vg + ^vq is /loo form bounded with relative bound 0. Recall that B > 
e implies a > 1. Choosing p = pc '■= —a 2 (^ + b) =: — a 2 v c will insure that 
(/ioo + ct 2 )rc(p)(hoo + a 2 ) 1 / 2 < 2 which is what we want to prove. A more careful 
argument, which we will skip, will yield an explicit b. QED 

We continue with the proof of Theorem 16.11 By l|113(l and i|112|) with p = pc, 
we find that 

if both 

a 3 / 2 o?l 2 
do® > ±C m v c ^ =: C' c ^ 

and a 1 / 2 ^ 1 / 4 > 4Cj72) which ' 

since B > a(B) 2 B is increasing, is equivalent to 
B > 4 3 C , 4 ^ q;(4C 2 ^ )~ 2 . Since we also need B > e we put 

43C fel 

Bp := max ^ — 2 , e 

4c/ 



This fixes our constants C' c and -B^, and also implies that ||r e ff(£)|| < 2||rc(£)|| 
by the resolvent formula. 
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To show that 1109fl implies that ||r"c(£)|| < 2||r e ir(£)||, we repeat the argument 
with rc and r c ff interchanged: by the resolvent formula, 

rc(0 = r cS (0 1/2 (l + KciO)' 1 kc ff (0| 1/2 , 

with 

K c (0 ■= \r eS (0\ 1/2 (h c - h cS )r cS (0 1/2 , 
so that, using lemmas [4.51 15.21 and !5. 31 we arrive at 

if > C^q; 3 / 2 ^^ 1 / 4 , by the first part. We therefore conclude that ||-Kc(£)|| < 

1/2, and hence \\r c (Q\\ < 2||r eff (f)||, if both ^B 1 ^ > 4Cjjm, which wiU bc 
satisfied if B > B' c defined above, and if 

%eff|CjJ2| 1 



y/ZBV*dc{&) ~ 2' 

The latter inequality is equivalent to 

/ N \ a 3 ' 2 

which yields condition H109(l . 
Finally, if £ satisfies then 

4C, |Za f "ca 2 . 

■=== — — — max < — — — — , 1 



||re ff (0-ra(OII < a i/ W *dc(0 ^ \d c (0 

a 3 ' 2 

a 3 ' 2 
C ' C d c (0 2 B^' 

where we used that <ic(£) < a 2 /2 < vqc? ■ This finishes the proof of Theorem 16. II 
QED 

Proof of Theorem M.cA We define 

(115) B c :=max{B eS ,B' c }, 
and 

(116) c c := max{^, 2c cS a(B c )- 1 / 2 B c lfi }- 
Suppose that B > Be, and that 

c c « 3/2 i?- 1/4 < d c (0 < \a 2 . 

By Theorem IO tZ e fr(£) < 2d c (£) < a 2 /2. By the same theorem, and by JTTSJl, 

4ff(£) > ido(0 > c^aiBcT^B-^aiB- 1 '* > c cS a{B)B- 1 ' 2 , 

since a(B)/^/B~ is a decreasing function of -B > 0. The conditions of Theorem ll.il 
are therefore met. Since the condition I|11U|) of Thcorcm l6.1l is clearly also satisfied, 
we conclude that the difference of the resolvents <|15|) can be estimated by 

a 2 , a 3 ' 2 n 3 / 2 

CcS- — 7= + C C , /Mini IA ^ ^ 
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for B> B c , with 

(117) C c :=±C eS ^^l + C > c , 

where we used that a(B)/y/~B is decreasing. QED 



7. Proof of Theorem 11.51 

This will be done by closely following the strategy of section 6. First, we compare 
the resolvents r e ff(£) of h e g and rg(£) ■= (hg — £) _1 of hg := hf' M . Recall, that 
d s (0 := dist (£,a(hg)). 

Theorem 7.1. Let a — a(B) be defined by Jj^). There exist (computable) constants 
Bg,C s > 0, only depending on Z, N and M, such that for all B > B' s and all real 
£ satisfying 

(118) d s (0 > C' s a, 

we have that £ € p{h c s), with ||r e ff(£)|| < 2||r<j(£)||. In addition, letting 
(H9) 

C s ' = mBx{C s , 4C(n2|(^ + l)} with C m := (nZ+ N{N ~ 1} ) C m , 

and £ e ff given by 1^96)) . then if 

(120) d s (0 > C'J a, 
we also have that \\rg(£)\\ < 2||r e ff(£)||. Finally, if 

(121) C' s a < d s (0 < \a 2 , 
then 

(122) \\r eS (0-rs(0\\<C' s " 



Proof. 

As in the proof of Theorem 15 . II 16 . ll one shows that 

ii (e\ < 1 11^(011 

\\ref£{Q - r s {Q\\ < 



ds(0 l-WKsiOr 

where Kg := Irg^)] 1 ' 2 ^^ — hg)rg(^) 1 / 2 . Using lemma l4~%l with c = 2, lemma 
14.61 the triangle inequality, and similar comparison arguments as in the proofs of 
Theorems I5.1l and l6.1l we easily find that, for any \i < inf a (kg), 



(123) uncoil < -GM max |_i|_,i| || r5(/i) i/ 2(ftoo + a 2 )1/2||2) 

where C jYT^l was defined above. We then use the following analogue of lemmas 15. 31 
andIO 

Lemma 7.2. Let vg = 1/2 + 4NZ 2 , fig := — vga 2 and let a be defined by Then 
fig < inf a(hg), and 

(124) ||r 5 (/i 5 ) 1 / 2 (/ l0 o + a 2 ) 1 / 2 || 2 <2. 
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Proof. We will use , as before, the scaling z i— > z/a. We get, with the help of lemma 

ED 

JV 

h$ — n ~ a 2 (/ioo — 2wa — oT 2 [i) > a 2 (hoQ — 2Z S(zj) — aT 2 n) 

> a 2 (/i 00 (l-2Ze)-7V£- 1 Z-a- 2 /i) (e > 0) 
= ^a 2 (h 0Q + l), 

if we choose e = e s := 1/(4Z) and ^ = /z 5 := -a 2 (^ + 47VZ 2 ). QED 

We continue with the proof of Theorem 17. II By 11123(1 and H124J1 with /i = /ij, 
we find that 

11^(011 <-^=* max |_ij<- 

if both 

d*(0 > 4Cjngj^ a =: C' s a, 
and a > 4C j^^[ which, since -B i— * a(-B) is increasing, is equivalent to 

B>B' S := IGC^e^m. 

This fixes our constants C' s and B' s , and also implies that ||7" e ff(£)ll — by 
the resolvent formula. To show that (|120|l implies that ||t\5(£)|| < 2||r e ff (£)||, we 
repeat the argument with r$ and r e ff interchanged: by the resolvent formula, 

rs(0 = r cS (0 1/2 (l + K s {0)~ 1 M0\ 1/2 , 

with 

K S (0 := \r eS (0\ 1/2 (h S ~ h eS )r eS (C) 1/2 , 
so that, using lemmas [4.81 15.21 and 15.31 we arrive at 

l|JWI1 - 2max i^)'7 — - 2max { w 1 / — ' 

if ds(£) > C^a, by the first part. We therefore conclude that ||-Ka(OII 1/2 if both 
a > 4C j^ i<j| , which is satisfied since B > B' s defined above, and if 

%eff|QTTTH < 1 

ad s (0 ~ 2' 
The latter inequality is equivalent to 

which yields condition (|1 20|> . Finally, if £ satisfies (|121|l , then 

where we used that ds(£,) < a 2 /2 < z^a 2 . This finishes the proof of Theorem 17. II 
QED 

Proof of Theorem \1.5\ We want to realize the assumptions of Theorem II .11 and 17. II 
with conditions on B and dg. We define 

B$ := maxjBcff,^}, 
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and 

(125) c s := max{C£, 2c cff F~ 1/2 }. 

Clearly B > B$ and cga < dg(£) < a 2 /A will do, for under these conditions, using 
Theorems Il.ll and l7.1l the left hand side of (|20|) can be estimated by 

CeffQ 2 C' 5 a ( AC cS a(B s ) A a = a 

since we know that both 2d c ff(C) > dg(£) x 

and 2ds(£) > d e s(£), by Theorem 17. II and since a(B)/^/B is a decreasing func- 
tion of B > 0. QED 



8. The fermionic case 

We first prove Theorem 11.61 This is simply done by repeating the proofs of 
Theorems 11.11 fOl and IT31 for the 'fermionized' operators, that is, for the operators 
sandwiched between P AS . We have to check that the main ingredients of these 
proofs remain valid. First of all, corollary 13. 31 is used to compare the resolvents of 
Hf' M and H^f + Wf. Next, lemma IP remains valid for W f := P AS WP AS , with 

the operator norm being the one on P AS ^Ran n^ M ^ , since P AS commutes with 

/ioo on Ran 1T^ M = L 2 (R N ,FS) (as we will explicitly see below, P AS not only 
mixes the coordinates of M. N , but also the different components with respect to the 
natural basis of F^: however, /loo acts in a scalar way). We then repeat the proof 
of Theorem 11.11 in section 5, replacing d c g everywhere by d e ff,f- Similar remarks 
apply to the proofs of Theorems 11.31 and 11.51 QED 

We next turn to Theorem 11.81 The parameter B here plays a non-essential 
role, and we will simply drop it, writing X m , Xm, F M for X^, ^m. f m> etc - We 
start by analyzing the subspace of anti-symmetric wave functions in the range of 
n Gff := Ilf ff . Recall that 

E(M) = {m = (mi, - • • ,m N ) : m,- > 0,mH \-m N = M}. 

The permutation group Sn acts on E(M) by ct-(toi, • • • , mjv) = (m 5 (i), • • • , m CT (An) 
and E(M) can therefore be written as a disjoint union of orbits of Sn '■ 

E(M) = (J S N • m, 

M C E(M) being a set of representatives of g \ E(M). If we let 

(126) Vfn = Span {X a . m : a £ S N }, 

then, recalling that Fm = Span {X m : to £ E(M)}, we have the orthogonal decom- 
position 

Fm = @ Fm- 

From this it follows that 

l4 M (L 2 (R 3N ) ® C 2N ) = L 2 (R N ) ® F M = L 2 (R N ) ® F w . 

Since P" 45 leaves each L 2 (R N ) ® F„ invariant, it suffices to analyze the subspace 
of anti-symmetric wave functions in each of the latter. We therefore fix an to £ M. 
and let 

(127) Gm — {a £ Sn : u -m — to}, 
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the stabilizer of to. Choose representatives a\, • • • , ok, K — A" (to), for the right 
equivalence classes of GW in Sn : Sn jGm = {fiGm, • • ■ , (JuGm) with <JiGm H 
(jjGm — if i 7^ j- Then X ai .m, ■ ■ ■ , X CK .m constitutes an orthonormal basis for 
Fjn, and each element ip = ip(x, y, z) of L 2 {R N ) ® can be uniquely written as: 

K K 

(128) ip = y^ajX^.jfj = y^ j aj(z)X aj . W i(x,y) ) 

3=1 i=i 

for suitable a 3 = a 3 (z) E L 2 (R N ). For any such ip and r = (a;, y, z) e M 3Ar = (M w ) 3 , 
i/j(t ■ r) = ^2 a j( r ' z)X aj .rn(T ■ x,t ■ y) 

3 

= ^2a,j{T ■ z)X( T -i aj) — L {x,y), 

3 

where t ■ r = (rym, • • • , r T (N)), and similarly for r • x, r • y and r • 2, and where we 
used that 

X m (r -x,r-y) = Y[xm j (x r (j),yrU)) = X r-^ m (x,y). 

3 

It follows that ip S L 2 (R N ) ® F„ is anti-symmetric iff, for any r S Sn, 
^a 3 (r • 0)X (T -i CTj) . m = (-l) T ^2a j (z)X IJj . m . 

3 3 

This is equivalent to the statement that ip is antisymmetric iff 

(129) ^a j {T-z)X {Ta . ) . m = {-iyJ2 a i( z ) X °i-™i 

3 3 

since the two statements are equivalent when r is a transposition, and these generate 
Sn- The version (|129fl . with no r _1 , will be more convenient to work with. We 
next observe that the map Uj — > raj gives rise to a a permutation p(r) of Sk '■ 

Lemma 8.1. For any to G M. the map p = Pm '■ Sn — > Sjf , A = K(jn), such that 

P( T )(«) = J ^ ro "i € CTjGW 
is a well defined homomorphism. 

In other words, p(r) is characterized by: 

rcr, e a p ( T ){i)Gm- 

One easily verifies that p(r) is indeed a permutation of {1, • • • , AT}, and that p is 
an homomorphism of Sn into S/^. 

With this notation, the left hand side of (|129fl reads : 

3 

and on replacing j by p(r)~ x {j) and using the fact that the X aj .jn form a basis of 
Fjn, we find that ip is anti-symmetric iff, for all j and all r € Sjy, 

(130) aj(z) = (-l) r a p(T -i )(i) (r • 2). 

If we successively replace r by r _1 and z by t • z, this becomes 

(131) Oj(r • 2) = (-l) T a p(T)U) (z), 

which implies that all a,j are uniquely determined by any one of them, a\, say, 
which we let, by definition, correspond to a\ = e, the unity element of Sn- More 
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explicitly, since p(<jj)(l) = j (for <Jj<J\ = <Jj € UjGwc = G p ( a . )(\)-m, by definition of 
p), equation l|130|) with r = <jj implies the important relation 

(132) a j (2) = (-l)% 1 (a r 4 

In particular, P AS (L 2 (R N ) ® f^) can be identified with a subspace of L 2 (R N ), 
by sending ip to ai (see ((13311 ) . We now analyze the symmetry properties of a\ 
imposed by the anti-symmetry of ip. 

Lemma 8.2. Let H be the subgroup of Sn generated by the set { Tcr ( ^)(i) ■ T ^ ^W}- 
Then, for all a € H , a\{a ■ z) = {— l) a ai(z), and these are the only symmetry- 
conditions which the anti-symmetry of ip imposes on a\. 

Proof. Let r G Sn be arbitrary. Then by 1(131(1 . ai(r • z) = {— l) T a p ( T )(i) (z) which, 
by JT33), equals (-l) T (-l) CT 'MWai(o>( T )(i) ■ z). Therefore 

whence the lemma. QED 



Lemma 8.3. The group H of lemma \H. 6 A is generated by the union of all stabilizers 
G ay w °f °~j - m > 1 <3 < K - 

Proof. Let r S Sn- Then r G OjGm, for some j. Since tcti = re G UjGm, we have 
that p(r)(l) = j, and therefore rcr^ )(1) G ajG^aJ 1 = G aj . m . 

Conversely, if a G G ay m, then er = crjer'cr,^ 1 , for some a' G Gm- Put r = cr/cr'. 
Then p(r)(l) = j, since o~^a'a\ G o-jGm, and therefore cr = rcr^ 1 = rcr p~( T )(i) is a 
generator of H. We conclude that the set of generators of H equals LljG ar rn, which 
proves the lemma. QED 

Lemma 8.4. Let H be the subgroup from lemma UTM Lf Gm — {e}, then H — {e}, 
while if Gm ^ {e} ; then H — Sn- 

Proof. It is obvious, from lemma that if G„ — {e}, then H = {e}, and a± does 
not have to satisfy any symmetry-conditions with respect to the action of Sn, by 
lemma 18.21 

Now SLippOSC tllclt GfrT 

is non-trivial. Then there exist two indices i and j such 
that rh~i = rh~j. We can suppose, without loss of generality, that i — 1 and j — 2. In 
that case, the transposition (12) is in Gm, and therefore (<r(l), er(2)) = a(12)a^ 1 G 
H, for all a G Sn, by lemma 1531 again. But then all transpositions will be in H, 
which clearly implies that H = Sn- QED 

Define a linear mapping 

(133) U m : P AS (L 2 (R N ) ® F w ) -> L 2 {R N ), 

by 

t / m('0)(2) = VK(^(-,z),Xrn) L 2 {R 2N, 

where we recall that K = K(m) = #{Sn /Gm). If ^> is given by 1(128(1 . then 
Um(ip)(z) = y/Ka 1 (z). By 1(131(1 . HC/wr^H 2 = ||V>|| 2 , so that J/™ is unitary and 
therefore injective. If Gm is non-trivial, then H = Sn, and the image of Um is con- 
tained in the space L\ S (M N ) of anti-symmetric wave- functions on M. N , by lemmas 
18.21 and 18.41 Since the only symmetries of a\ are those imposed by H, it follows 
that Um is surjective onto L\ S (R N ). Similarly, if Gm = {e}, then H = {e}, and 
the image of Um is L 2 (M. N ): indeed, if a\ = a\(z) is arbitrary, then 

ip ai := ^2 ai(a- z)X a .rn(y,z) 

cr£S N 

is an anti-symmetric element of L 2 (M. N ) gj Fm such that Um(ip ai ) = oi. 
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Proof of Theorem ] 1.81 Recall that 

Mi = {meM:G 5f ={e}} I M 2 = {m e M : Gjn £ {e} }, 

and define 

e e 

(134) U M =:U^:P AS {L 2 (R N ®F M )) -> ]T L 2 (R N ) © £ L^fM"), 

by /7m := $seM^m- Then we have shown that ?7m is a surjective isometry. The 
intertwining formula of h^ { with Um being obvious, this proves Theorem 1 1.81 QED 

Contrary to hg { , the operator fM^cf^M wn ^ m general not act diagonally any- 
more on the range of Um, but will contain terms which couple anti-symmetric and 
boltzonic components in (|134fl . that is, components in L^ s (R Ar ) and L 2 (R N ). The 
potentially problematic terms in h^ f are UmC^^U^ and UmP^^U m (still drop- 
ping the B from our notations). The first one is easily seen to act diagonally on 
the right hand side of (|134f> : recall that 

f~tn,M _ 



and identify this with a Hermitian operator on Fm — Span {X m : m G S(M} (it 
acts as a multiplication operator on L 2 (R N , Fm)). The matrix of C"^f in the basis 
X m , m € £(M), is easily seen to be diagonal. Moreover, 

(X ff . m |0*™) = -/ \X m (a- 1 ■ (x,y)\ 2 log\ {nf =1 p 2 ) 1/N dxdy 

= (X m |C™^|X m ). 



Hence, taking B — 1, for simplicity, C™^ will act on f„ as scalar multiplication 
by 



( Xw p:^\x m ) = (log 4) nf =1 i4i 2 (p„)^4 



1 N 

-Y 



-iE( 2 ^ ! ) / P 2 ^ +1 ^g4e- p2/2 dp 



3 = 1 



1 



JV 




log2 ~ it \ H^ m 

where ip(z) = T' (z)/T(z), the logarithmic derivative of the T-function. Hence 
C%C'lv i^M simply acts diagonally on the right hand side of ^ 134(1 . and more pre- 
cisely, by scalar multiplication in each component. 

The operator UmC%^U m is more complicated: it is not going to be diagonal 
in the natural basis and it will in general mix the different components Ran Um, 
even those with index in A4i and M.2- This is already the case in the simplest case 
in which both M.\ and M.2 are non-empty, namely that of two electrons, N = 2, 
and a total angular momentum of M = 2. In that case £(M) is the union of two 
orbits under S2, namely {(0, 2), (2, 0)} and {(1, 1)} the first having as stabilizer the 
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identity, and the second having as stabilizer the full group 52- We can therefore 
take Mi = {(0, 2)} and Ai\ = {(1, 1)}. We will now compute the matrix element 

( l35 ) (^(0,2) 1^:21^(1,1))' 

and simply observe that the result is non-zero. For this computation it is convenient 
to use complex notation for the lowest Landau functions (j25(l : if we let C = % + iy, 
then (taking B — 1 again) 

y 1 - c r i e - |c|2/4 

Am ^rn^ c , 

where c m = (2ir2 m m\)~ 1 / 2 . We then find that l|135|) equals 



Making the (by now familiar) change of variables u = ((1+^2) /V2, v — (Ci — C2)/v / 2, 
this integral becomes 



c c 2 c 2 



4 



(u 2 -v 2 )(u~v) 2 log(M 2 /8)) e-d u l +1*1 ^ 2 dudv. 

IC JC 

Since (u 2 — v 2 )(u 2 — 2uU + v 2 ) = \u\ 4 — \v\ 4 + 2i Im v?v 2 + 2Wu(v 2 — u 2 ), and since, 
in general, 

/ [ u a u (i v u v K g(\u\ 2 ,\v\ 2 )dudv = 0, 
Jc Jc 

unless a = (3 and v — k, we obtain that (|135|1 equals 



cpc 2 c\ 

„2 



{\u\ A - \v\ A )\og{\v\ 2 l$)e-^+^l 4 dudv 



C JC 



.£°^£l / | u | 4 e-H 2 /2 du / l g(| w | 2 /8) e -H 2 / 2 ^ 
4 [Jc Jc 

e-^ /2 du [ |w| 4 log(|w| 2 /8)e- |l,|2/2 d« 



16^2 

as can be shown using the T-function and its derivative. So (|135|) is non-zero, and 
C^M^cf^M mrxes the two sectors. 

9. Concluding remarks 

We finally want to give an idea what these effective Hamiltonians might be good 
for. Our original motivation for introducing them was for studying the structure 
of the bottom of spectrum of H B ' M , in connection with the maximum ionization 
problem for atoms in strong magnetic fields; see below. To illustrate how this works, 
we first consider the comparison withfff M = hf' M ®H» < M ; cf. TheoremO bmce, 
by construction, B$ > -Sj^gJ, we know from Theorem 13 . II that a (^(H^ M ^j C (0,oo) 

if B > Bg. Let Eg := inf /if, then clearly Eg < 0. 

Since hf is unitarily equivalent to a 2 (—^A z +2vg) we conclude that whenever Eg 
is an isolated eigenvalue, its position as well as its isolation distance is proportional 
to a 2 . That Eg is an eigenvalue is true for Z large enough when N is fixed; to 
determine how big Z has to be exactly, relative to N, for this to happen is an open 
problem (see below). Let us assume henceforth that E$ is an eigenvalue, which then 
necessarily is simple. Consequently Eg is an eigenvalue of Hg' M with multiplicity 
dimF^, see JTHJl. Choose two points £± := E s ± Co? in p{Hf M ) n p(H B < M ), 
which satisfy H19|) . This is possible when B > Bg, see Theorem 11.51 Let L be 
the circle in the complex plane centered at Es with radius Co? and define P and 
Ps as the eigenprojections associated to H B ' M and h B ' M , respectively, onto their 
spectrum inside F. To estimate P — Pg we need a bound on i?(£) — Rs(0 ■— 
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(tf^-O^-Cfff" -O -1 for all ^ G T. We know already that \\R(Z±)- Rs{0\\ < 
Csa/ds(£,±) 2 = CsC~ 2 a~ 3 by theorcm ll.5l To propagate this estimate on all of T 
we use the convenient formula (see [K] IV. (3. 10)]) 



\R(0-Rs(0\\< 





B,M 

H s ~ 5 


2 

ll««±)-fl«(C±)ll 


i-|f-£±l 




ll««±)-«s(f)ll 



< 



C s C-- 1 a- i 



Then integrating over the contour T finally gives \\P— Pg\\ = 0(a~ 1 ) as B tends to 
infinity. This shows that for B large enough these two projections have the same 
dimension and since they are continuous with respect to B we finally get that for 
all B > B s , dimP = dim Pa = dimP^. 

Our conclusion is therefore that, for sufficiently large P, H B ' M will have a cluster 
of eigenvalues in the interval (Es—csa 2 , Eg+cga 2 ) with total multiplicity of dim F^, 
and apart from this no eigenvalues at a distance Ca 2 from E$, (the allowed P's 
will depend on C), so that the cluster is separated from the rest of the spectrum 
of by a distance proportional to a 2 . In the particular case when dimP^ = 1, 
i.e. M = or N = 1, we get an estimate on the difference of the eigenvectors 
<f> — $5 = 0(a^ 1 ) as B — » oo ( in the L 2 norm). Here <& denotes the eigenvector of 
h b,m and $ s ( Xt y } z j ■— i fi f( z )X B =0 (x, y) if M = 0, where ipf is the ground state 
of hf. In the case N = 1 one has tpf(z) = \j2aZ e ~ 2aZ ^ . 

Consider now the comparison with h^ ,m . The difficulty here is to find the nec- 
essary a priori information on the structure of a(h^' M ). In the case of N = 1, 
using the invariance of a(h^' M ) under the reflexion z i— > —z, and the characterisa- 
tion of the domain of <r(h^ M ) in Appendix A, one sees that the odd spectrum 6 is 
P-independent and coincides with the spectrum of the hydrogen in the s sector of 
symmetry. Since the even spectrum intertwines with the odd spectrum and since 
it is monotonically decreasing with respect to P, cf. (|ll|l . it is easy to realize that 
cr(h c ' ) l~l K_, apart from the ground state energy, is made up of clusters of two 
eigenvalues, the clusters being separated by a distance of order 1 as P — ► oo. Thus 
by using Theorem 11.31 and following a similar strategy as above one can conclude 
that for N = 1, and arbitrarily small e, a(H B ' U ) n (— oo, — e] has the same cluster 
structure for P > P e sufficiently large. Moreover this spectrum deviates from the 
one of the Coulomb model by at most cca(B) 5 P~4 as P — * oo. 

The model operator /i^ M is for the moment of mainly theoretical interest since 

it does not seem to be solvable even in the one electron case. Notice however 

that one could solve /i^ M numerically, at least for few electrons and small M, 

and subsequently use Theorem 1 1.1 1 to approximate the true spectrum of H B ' M for 

large P. Given the non-trivial dimension reduction achieved by theorem 11.11 ffrom 

wave-functions of 3iV variables to ones of variables, albeit vector-valued) such 

a procedure would, from a numerical point of view, seem preferable to attacking 
H B,M directly . 

Whether the simpler models, i.e. the delta and the Coulomb model, are solvable 
in the A-electron case, N > 2, is a challenging question in view of applications. 
We are thinking in particular of the problem of determining the maximum number 
N c of electrons which a clamped nucleus with charge Z can bind when an intense 
homogeneous magnetic held is applied. LSY has shown that liminf N c /Z > 2 as 
Z, B/Z 3 — > oo. Very little precise is known for fixed Z and high P. It is conjectured 
that there should be a P-independent absolute (that is, non- asymptotic) upper 



6 that is, the spectrum of he restricted to the odd wave functions 



EFFECTIVE HAMILTONIANS FOR ATOMS IN VERY STRONG MAGNETIC FIELDS 39 



bound of the form N c < aZ + b, similar to Lieb's bound N c < 2Z + 1 valid when 
-6 = 0, but this is as yet unproved. Some weaker results are known, of which the 
best to date is the one of |Sei| : see also |BR| for work on heuristic models related to 
our h A . It is natural to first try to solve the maximal binding question for hf , or 
any of our other effective Hamiltonian, and use the approximation theorems of this 
paper to draw conclusions for H B ' M itself. Some modest progress is possible in this 
way. It is for example known that the delta model with two electrons is at least 
numerically solvable, see |Ros| . and that this model possesses a unique bound state 
at the bottom of its spectrum as long as Z > 0.375. Therefore using Theorem 1 1.51 
we see that for all Z > 0.375 there exists Bz > such that for all B > Bz one 
nucleus with such a charge can bind two electrons. As a consequence, Lieb's bound 
of N c < 2Z + 1 is no longer valid in strong magnetic fields. For general Z, no 
maximum ionization bound for the (5-model is known as yet. 

Let us now briefly turn to the effect of particle symmetry. It follows from The- 
orems II. 61 and 11.81 that H^' M can be approximated by a direct sum of copies of 
hf acting on anti-symmetric L 2 (R N ) plus a direct sum of copies of hf acting on 
L 2 (R N ) without any symmetry condition. The latter will occur iff A4i ^ 0, which 

is the case iff M > 0+H h(iV-l) = ^N(N-l). For such M, the ground state 

energy will be approximately that of boltzonic hf , which is also the ground state 
energy of the bosonic hf , and the same can be shown to be the case for the ground 
state wave function (assuming there is one), by standard permutation arguments. 
In fact, Theorem 11.61 plus Theorem II .81 predict the existence of a cluster of #Aii 
eigenvalues at the bottom of the spectrum of H^' M , at a distance of order a(B) 2 
from the origin as B — > oo. 

A further interesting corollary to Theorem 11.81 can be obtained by considering 
M as a free parameter. If hf possesses a ground state, whose energy is isolated in 
its spectrum, then for sufficiently large B, Hf := p AS H B will assume its ground 
state for an M > ^N(N — 1). Stated otherwise, assuming there is a mechanism 
for transfer of the angular momentum (e.g. emission and absorption of photons), 
atoms in strong magnetic fields will have an orbital angular momentum in the field 
direction of at least ^N(N — 1). A natural conjecture is that we have equality here. 
Notice that this conjecture was shown to be true in the case of N = 1 in AHS , 
see also jBaSe) . 

We mention one further application of these effective Hamiltonians. After the 
location of the spectrum to leading order one can now use regular perturbation 
theory to compute lower order corrections. We have shown how this can be done 
m BcBDP . This seems definitely more convenient than variational techniques 
and more familiar than the Birman-Schwinger method used in |AHS| for the one 
electron case. Continuing for example the above comparison of H B,U with h s ' , it 
is immediate to realize that adding the first order perturbative correction will give 
an error of order 1. In case N — 1 we get that the ground state energy of H B,M 
is equal to (hff 1 ^,^) + 0(1) with $ s (x,y, z) = V2c7Ze- 2aZ W x&(x,y). This 
should be compared to Theorem 2.5 of |AHSj . In fact, one can write the ground 
state energy as a convergent power series, each term of which being of order a~ k 
for k running from —2 to infinity. However, as pointed out in [AHS , in view of the 
log(i?) behaviour of a this series is of limited value. The situation will be much 
better with the Coulomb model h^' M since the perturbation series will converge 
much faster because of the a^B~i behaviour of the r.h.s. of l)15|l. We hope to 
tackle this program soon. 

The above remarks concentrated on applications of our effective Hamiltonians 
to the ground state energy of H B,M , but their potential interest is not limited to 
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that. As the example of he shows, other parts of the discrete spectrum of H B,M 
will also become amenable to analysis, if this is the case for the effective operator. 
On a conceptual level, Theorem 11.11 gives a precise mathematical sense to, and 
justification of, the physicist's attractive heuristic picture of an atom in a strong 
homogeneous magnetic field as consisting of electrons in their lowest Landau band 
states interacting through a kind of " residual" electrostatic interaction. Finally let 
us note that the technics developped in this article are expected to work in other 
contexts. An interesting example is that of 2-dimensional electronic systems on a 
cylinder which describe excitons in carbon nanotubes; cf. [CDP . 



Appendix A. A characterization of the operator domain of he 

In this appendix we will characterize the operator domain of he = h B . It will in 
fact be convenient to consider a slightly more general situation. Let C — {£„, 1 < 
v < K}, be a finite collection of hyperplanes in R N (we might more generally con- 
sider non-singular C^-hypersurfaces). Let F be a finite dimensional complex vector 
space, with Hermitian inner product (•, •) and let A v , B v be Hermitian operators on 
F. Let H k (R N ,F) be the fc-th Sobolev space on R N , with values in F. Then the 
following sesquilinear form is well-defined on H (R , F) and bounded from below: 

(136) t c (u,v) = i(Vu, Vv) + ]T <L*Pf (| • r 1 ), (A u u,v)) + (L*J, {B v u,v)) , 
(•, •) denoting the duality between distributions and test functions. We let 

he = -\& z + £ A v Pf + B v 5{L v (z)). 

be the associated self-adjoint operator, whose existence is guaranteed by the KL 2 MN- 
Theorem. Let TZ be the set of connected components of R^ \ U„Ker L V: so that 
\ U„Ker L v = U Re uR and R n R' = 0, if R, R' G TZ, R ^ R' . Note that, on any 
of these components R, L*Pf(l/| • |) simply equals the function l/\Lj(z)\. We will 
identify Lj with an element of R , using the Euclidean inner product on R^; the 
latter will be denoted by a dot: z • w, to distinguish it from the Hermitian inner 
product (v,w) on F. If we let H 2 {R 1 F) be the F- valued Sobolev space of order 2 
on the open subset R C R^, then the domain of he can be characterized as follows: 

Theorem A.l. Let u G L 2 (R N ,F). Then u G Dom(hc) iff the following three 
conditions hold: 

(i) u G H l (R N ,F), 

(ii) For each R G TZ, 

(in) For each j and for each x G {Lj = 0} \ Uk^j{L)~ = 0}: 

Lj ■ Vu (x + e J^j^j - L j • Vu M5 - e jf^J — -HogeAju(x) + 2Bju(x) + o(l), 
as e — > 0. 

Proof. We will use the following characterization ( see |K]) of Dom(hc)' 

u G Dom(t £ ) = H 1 (R N ,F), 

(137) u G Dom(/i £ ) { and 

\tc(u,v)\<C u \\v\\ 2 , V W GDom(t £ )=iJ 1 (R Ar ,F), 
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(139) (L*Pf (1/| ■\)M = -I a51 ^T/ l 2 BI V " L ■ Vu(z)dz. 



the norm on the right being the L 2 -norm. Here we may, and will, suppose without 
loss of generality that v G C\ (M N , F) , the space of compactly supported F-valucd 
C 1 -functions. 

To analyze Ijl37|) . we will first establish a convenient integral expression for the 
pull-backs of 5 and of Pf (1/| • |) under a linear map L : R N — > M. Using the 
Euclidean inner product, we can identify L with an element of M. N , which we also 
denote by L. Using the definition of pull-back (cf. Hormander jHoj . chapter 6), 
one easily shows that if do^ denotes the Euclidean surface measure on Ker(L), and 
||L|| is the Euclidean norm of L, then 

(138) (L*5, ip) = f jr^da L , 

J{L=0} \\ L \\ 

and 

S ga(L(z)))log\L(z)\ 
\L\\ 

In fact, since taking pull-backs is coordinate invariant, it suffices to verify these 
formulas in a orthogonal coordinate system in which L — (\\L\\, 0, • • • ,0). 

We next observe that, for each R G H, there exists a (unique) function sr : 
{1, • • • , K} -> {0, 1}, such that 

R = {z G R N : (-ly^Ljiz) > 0,j = 1, ••• ,K}, 

and if e > 0, we define R e by: 

Re = {zeR N : {-ly^Ljiz) >e,j = l,---,K}, 

Observe that the boundary of R £ , as well as that of R, are polyhedrae, each of 
whose faces are contained in one the hypersurfaces {(— l) SR ^'Lj — e} and {Lj = 0}, 
respectively. On such a face, the outward-pointing normal nn. £ can be identified 
with the vector — (— 1) SR ^ Lj/\\Lj\\ (translated to the relevant base point on the 
face, to be precise). Note, that n^, e is only defined a.e. on the boundary (with 
respect to the surface measure). This will not cause difficulties, though. 

Suppose now that u G Dom(hc) and in particular satisfies the conditions i|137fl . 
We then have, for any v G C\ (R N , F) , 

(140) t c (u,v) = limJ^^ f {Vu,Vv)dz 
£ ^ { Ren ■'^ 

- 2^}^ rn- j|2 {L v -W)(A v u,v) dz 

- R Jr. Il-M 

^J{L„=0} \\l>v\\ J 

We want to apply Gauss' divergence Theorem to each of the integrals over R e : this 
is allowed since if u G 7? 1 (M Ar ,F) satisfies l|137|> . then by choosing v compactly 
supported in R, we see that 



2 * ^ " \L(z) 

which obviously implies that Au G L^ oc (R,F) (since we are staying away from the 
singularities on the hyperplanes), and hence u G H^ oc (R,F). Now 

(141) / (Vit, Vv)dz = I (-Au,v)dz + / (n RtS ■ Vu,v)da R . e , 

Jr. Jr. Jor. 
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dan j£ being the euclidean surface measure on the boundary. Furthermore, 
-||L,||- 2 / flE Sgn (L u ) log \L„\ (L u -V)[(A v u,v)} dz = 
f (A u u v) f 

JR E \ L v\ z )\ JdR e 

Using this, the second term in l|14(Jf> can, after re-arranging, be written as: 

(Aju, v) 



Ren 



+ EE ±H^II _1 f Lj=±s , sgn(L,(z)) log|i,(z)| (A jU ,v) 

j ± \l,A > e'it v ? j 

+ E ill^H" 1 L h =± e , s MLj(z)) log\Lj(z)\ (A jU ,v), 

j,k:j^k ± l-M > e if v 4. k 

the surface integrals being with respect to the natural surface measures. Since v 
is compactly supported, the third sum will vanish, in the limit of e —> 0. This 
follows from the local integrability of u ■ log | ij | on {L^ = 0} (if j 7^ k), which 
can be seen as follows. Since u £ H 1 , its restriction u|{£ fc= o} is in (vector-valued) 
L 2 . On the other hand, log|Lj|, restricted to = 0} is in L 1 2 oc (R Ar_1 ), and 
therefore itlog|Lj|, restricted to {Lk = 0} is locally integrable. Rearranging also 

the sum over 1Z of the boundary terms in (|141|) as a sum of integrals over the various 
hypersurfaces {Lj = ±e}, we conclude that for v £ C^.(R N , F), 



■ f ( 1 A V- A » U \ 

im / Au + > — — - ,v 

-^oJ Re [ 2 ^\L v {z)\ ) 



tc(u,v) = lim / ( -^Au + ^ rr7 -^ , v j dz 

+ E^E/ %nfW^> 

+ / II T II dcr {^=0}- 

J{L,=t)} \\ L j\\ 

Since u £ H 1 implies that + EZy/||Zy||) — u(-)IU 2 ({l =0}) = 0(y/e), and since 
(u, v) £ L ({Lj — 0}), we can replace the next to last term by 

21oge|| J L J |r 1 / (A jU ,v)da {L 0} + o(l), s -» 0. 

•/{Z. 3 =0} 

We can, for the same reason, replace v in the second integral by its restriction to 
{Lj = 0}. Letting e — > 0, it follows that if it satisfies Ijl37|l . then it satisfies (i), (ii) 
and (iii). 

Conversely, suppose that u satisfies conditions (i), (ii) and (iii) of the theorem. 
Then running the argument backwards shows that \tc(u, v)\ < C u ||w|| 2 , first for all 
compactly supported v and thence for all v £ 7? 1 (M Ar , F). Hence u £ Dom(lic), by 
(HSZI- QED 
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